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PPSFACE 


As  part  of  the  ILS  Performance  Prediction  program  (PPA  No. 
FA307) ,  a  first  phase  ILS  Localizer  performance  prediction  computer 
program  package  has  been  prepared.  This  package  consists  of  the 
computer  program  and  the  present  document  which  describes  the 
capabilities  and  limitations  of  the  computer  model  as  well  as 
the  step  by  step  running  of  the  computer  program. 

The  computer  program  is  intended  as  an  aid  in  predicting  the 
performance  of  different  ILS  Localizer  antenna  candidates  for  a 
proposed  runway  instrumentation  or  for  the  upgrading  of  an  already 
instrumented  runway.  It  is  also  intended  to  provide  a  relatively 
inexpensive  means  by  which  the  effect  of  any  proposed  changes  to 
an  airport  environment  (addition  of  terminal  buildings,  hangars, 
etc.)  on  ILS  performance  may  be  predicted. 

This  document  was  prepared  for  TSC  by  D.  Newsom  assigned 
full  time  as  a  programmer  to  the  ILS  Performance  Prediction 
program  and  by  A.  Watson  who  helped  in  its  writing.  The  document 
and  attached  computer  program  are  based  on  the  theories  and 
analyses  developed  by  the  TSC  group  (Chin,  Jordan,  Kahn  and  Morin) 
for  the  ILS  program  sponsored  by  H.  Butts  of  the  Systems  Research 
and  Development  Service  of  the  FAA. 
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1.  DEFINITION  OF  INSTRUMENT  LANDING  SYSTEM 


The  ILSLOC  program  has  been  written  to  simulate  certain  air¬ 
port  conditions  which  affect  the  localizer  portion  of  the  Instrument 
Landing  System.  The  ILS  is  used  to  provide  signals  for  the  safe 
navigation  of  landing  aircraft  during  periods  of  low  cloud  cover 
and  other  conditions  of  restricted  visual  range.  Separate  systems 
are  used  to  communicate  vertical  and  horizontal  information;  the 
horizontal  system  is  called  the  "localizer”. 

This  system  operates  by  the  transmission  of  an  RF  carrier, 
amplitude  modulated  by  two  audio  frequencies,  beamed  to  approaching 
airborne  receivers.  In  an  instrumented  aircraft,  the  localizer 
receiver  serves  to  demodulate  the  RF  signal ,  amplify  and  isolate 
the  corresponding  audio  signals  and  derive  a  signal  to  drive  the 
ILS  horizontal  display  in  the  cockpit.  The  pilot,  by  reading  the 
display,  can  determine  if  he  is  on  course,  to  the  left  of  the 
runway,  or  the  right  of  the  runway.  These  signals  must  be  strong 
enough  to  cover  a  radius  of  twenty-five  miles  around  the  antenna. 

The  directional  information  is  determined  by  the  relative 
strengths  of  the  transmitted  sideband  signals.  The  audio  frequency 
modulations,  which  are  fixed  at  90  H  and  150  Hz,  are  radiated 
in  different  angular  patterns  with  respect  to  the  runway  centerline 
extended.  The  "course"  is  defined  as  the  locus  of  points  where  the 
amplitudes  of  the  two  modulations  are  equal.  The  display  of  a 
difference  of  the  amplitudes  (90  Hz  and  150  Hz)  of  the  sidebands 
is  referred  to  as  the  Course  Deviation  Indication.  Thus,  the 
CDI  is  the  pilot's  indication  as  to  what  his  bearing  is  relative 
to  the  center  line  of  the  runway.  The  CDI  is  measured  in  microamps. 
The  actual  course  generated  by  any  particular  ILS  installation  will 
deviate  from  the  ideal  due  to  the  interference  of  spurious  re¬ 
flections  from  buildings  present  in  the  range  of  the  transmitting 
antenna.  The  deviation,  caused  by  these  buildings,  or  scatterrrs 
of  the  CDI  from  what  the  receiver  should  read  ideally  at  that 
point  in  space  (e.g.,  on  the  center  of  the  runway  and  CDI  reading 
other  than  0)  is  the  derogation  effect. 
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The  Localizer  system  transmits  an  asymmetrical  pattern  by 
beaming  a  "carrier  plus  sideband"  pattern  and  a  "sideband  only" 
pattern,  the  composite  of  which  gives  the  desired  effect.  If  a 
specific  localizer  system  uses  two  antenna  arrays,  four  sets  of 
signals  will  be  transmitted;  if  the  system  uses  a  single  antenna 
array,  two  sets  will  be  transmitted. 
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2.  ANTENNA  PATTERNS 


The  proper  angular  variation  of  the  transmitted  90  Hz  and 
the  150  Hz  modulation  is  achieved  by  the  radiation  of  two  independent 
sideband  patterns  by  the  transmitting  antenna  arrays.  Equal 
magnitudes  of  90  Hz  and  150  Hz  modulation  are  transmitted  in  each 
of  these  patterns,  however  with  different  relative  phases.  One 
of  the  patterns  is  symmetrical  with  respect  to  the  prescribed 
course.  An  unmodulated  carrier  wave  is  transmitted  with  the  same 
pattern  and  the  combination  is  commonly  referred  to  as  the  "car¬ 
rier  plus  sidebands"  (C  ♦  S)  signal.  The  other  signal  is  trans¬ 
mitted  in  an  "anti -symmetrical"  pattern  and  is  referred  to  as  the 
"sidebands -only"  signal. 

Figure  1  illustrates  how  these  features  are  used  to  obtain 
the  desired  directional  CDI.  The  magnitudes  of  the  C  +  S  and  SO 
sideband  patterns  as  functions  of  angu  nr  deviation  from  the 
course  are  illustrated  in  Figures  la.  The  sideband  amplitude  of 
the  C  +  S  pattern  represents  201  modulation  of  the  carrier  wave 
(or  a  "depth  of  modulation"  of  0.2)  at  both  SO  Hz  and  150  Hz. 
Considering  the  phases  of  both  modulations  of  the  C  +  S  signal  to 
be  positive,  the  relative  phases  and  typical  amplitudes  of  the  two 
SO  modulations  are  as  shown  in  Figures  lb.  The  resultant  90  Hz 
and  150  Hz  modulation  patterns  in  the  total  ILS  signal  are  obtained 
by  algebraically  combining  the  respective  C  +  S  and  SO  sideband 
patterns  (Figures  lc).  The  evident  consequence  is  that  the  depth 
of  modulation  is  greater  for  90Hz  than  for  150  Hz  to  the  left 
of  the  course  as  seen  from  an  approaching  aircraft,  and  the  op¬ 
posite  is  true  to  the  right  of  the  course.  This  difference  when 
properly  calibrated  in  relation  to  the  total  modulation  (90  Hz 
+  150  Hz)  reaching  the  aircraft  rfceiver  gives  the  CDI  as  appears 
in  Figure  Id. 

Since  the  strength  of  C  *  S  and  50  signals  fall  off  at  the 
same  rate  with  distance  from  the  transmitting  antenna,  the  CDI 
is  independent  of  range. 
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Figure  lb  Relative  Amplitudes  and  Phases  In  SO  Pattern 


Figure  lc  Resultant  Modulation  Patterns 


Figure  Id  Course  Deviation  Indication  (CDI) 
Figure  1.  Antenna  Patterns  Sketch 


FAA  standards  for  the  ILS  specify  that  within  a  certain  nar¬ 
row  angular  range  about  the  course,  vhe  CDI  should  be  closely 
proportional  to  the  aircraft’s  angular  deviation  from  course.  This 
sector  near  the  ideal  approach  is  termed  the  "course  sector"  and 
usually  extends  between  1  1/2°  and  3°  to  either  side  of  the  runway 
centerline.  The  wider  sectors  on  either  side  of  the  course  sector 
are  called  the  "clearance  sectors".  In  these  sectors,  which  extend 
a  minimum  of  35°  from  the  course,  the  CDI  is  required  to  always 
exceed  a  certain  minimum  magnitude.  The  presence  of  structures 
in  the  clearance  sectors  which  scatter  spurious  signals  into  the 
course  sector  is  the  primary  cause  of  derogation  of  the  localizer 
CDI.  Such  structures  are  illuminated  by  carrier  and  sideband 
signals.  The  ratios  of  150  Hz  modulation  to  90  Hz  modulation  in 
these  signals  are  determined  by  the  angular  position  of  the 
structure  with  respect  to  the  runway.  In  general  these  ratios  are 
different  from  those  transmitted  toward  the  aircraft,  due  to  the 
difference  in  angular  position.  The  signals  transmitted  toward 
the  scatterer  will  be  reflected  toward  the  aircraft.  Thus  the 
aircraft  will  receive  the  summations  of  the  direct  and  scattered 
signals.  Since,  in  general,  the  scattered  signals  will  have  im¬ 
proper  ratios  their  effect  is  to  distort  the  CDI.  To  combat  this 
problem  several  new  antenna  systems  have  been  designed.  Two  basic 
systems  are  used:  the  single  antenna,  and  the  "capture  effect 
system." 

The  single  antenna  system  radiates  two  patterns  from  one 
antenna  array.  The  signal  generated  in  the  course  sector  is 
stronger  than  that  generated  in  the  c7 sarance  sector.  However, 
because  of  the  derogation  effects,  the  signals  are  often  not  ac¬ 
curate  enough  to  meet  category  II  or  III  requirements  and  the 
more  accurate  "capture  effect  system"  is  used.  This  system  uses 
one  antenna  array  to  broadcast  a  very  narrow,  powerful  beam  in 
the  course  sector.  The  second  antenna  array  broadcasts  a  broader 
pattern,  at  a  slightly  different  carrier  frequency,  which  covers 
the  clearance  area.  This  system  diminishes  the  derogation  effects 
because  of  the  dual  frequency.  The  terra  "capture  effect"  has 
been  used  to  describe  this  two  antenna  array  system  because  the 
airplane  receiver  is  "captured"  by  the  stronger  transmission  signal. 
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3.  ILS  SIMULATION  DESCRIPTION 


The  ILS  simulation  program  makes  it  possible  for  airport 
planners  to  determine  what  the  effects  of  potential  airport 
buildings  on  the  ILS  performance  are  going  to  be.  Thus,  for 
example,  if  a  new  terminal  or  hotel  is  planned,  the  information 
as  to  size  and  location  of  the  building  can  be  input  to  the  program 
and  the  derogation  effect  of  that  building  can  be  determined. 
Because  the  derogation  effect  of  these  scatterers  is  so  important, 
the  program  can  warn  the  planner  ahead  of  time  to  change  the 
orientation  or  location  of  the  building,  or  it  can  assure  him 
that  the  building  would  not  jeopardize  the  airport's  current 
FAA  rating. 

The  output  of  this  program  is  a  magnetic  tape  of  values  of 
the  CDI.  Graphs  are  generated  by  a  plotting  routine  (using  the 
values  derived  from  the  ILSLOC  program)  to  show  the  CDI  in  micro¬ 
amperes,  along  a  flight  path,  for  the  scattering  surfaces  input. 
These  generated  graphs  would  serve  the  same  purpose  as  the  FAA 
strip  charts  which  are  generated  for  a  certifying  flight.  The 
simulation  graph  differs  from  the  actual  recorded  measurements 
due  to  limitations  of  the  program  which  will  be  explained  later 
in  the  text. 

The  ILSLOC  program  simulates:  transmission  from  the  various 
types  of  localizer  antenna  systems;  the  trajectory  of  an  aircraft 
flight  over  which  the  CDI  is  to  be  determined;  and  the  scattering 
from  rectangular  and  cylindrical  surfaces.  The  program  permits 
various  simulated  flight  paths. 

The  program  is  not  an  exact  simulation  of  the  certifying 
flight,  due  to  certain  simplifying  assumptions  which  were  made. 
These  assumptions  include: 

a.  A  flat  perfectly  conducting  ground  plane 

b.  Perfectly  conducting  reflectors 
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c.  Far  field  scattering  -  ali  scattering  from  a  surface 
is  assumed  independent  of  all  other  surfaces,  thus 
multiple  reflections  from  walls  and  near  field 
interactions  are  ignored. 

d.  A  noise  free  environment 

e.  Relative  field  strengths  -  the  absolute  field  strengths 
involved  are  not  calculated.  Thus  while  we  can  calculate 
the  CDI's  in  microamperes  we  do  not  ascertain  the 
absolute  electric  field  intensities. 

f.  An  idealized  1LS  receiver  model. 

In  addition  to  these  assumptions  the  approximations  of  the 
scatterer  can  lose  accuracy  when  the  dimensions  approach  less 
than  a  few  wavelengths.  Since  the  program  determines  the  scat¬ 
tering  from  a  surface  independently  from  all  other  scatterers,  the 
shadowing  of  one  structure  on  another  is  not  included.  Thus  if 
one  building  is  between  the  antenna  system  and  another  building, 
it  will  shield  the  second  one  from  some  or  all  of  the  ILS  signal. 
The  amount  of  energy  reaching  the  second  building  will  depend  upon 
diffraction  effects  which  are,  in  general,  too  complicated  to 
analyze.  It  may  be  noted,  however,  that  diffraction  effects 
themselves  are  included  as  part  of  the  physical  optics  approxima¬ 
tion  used  (Ref.  1).  By  using  rule  of  thumb  approximations  the 
analyst  can  determine  roughly  how  much  power  will  reach  the 
second  building.  If  the  level  is  small  the  building  may  be 
ignored  completely.  If  on  the  other  hand  the  power  level  is 
large  then  the  structure  should  probably  be  included  as  though 
there  was  no  shielding  effect.  This  will  give  a  conservative 
CDI  estimate  (i.e.  larger  derogation  than  actual),  but  this 
will  serve  for  most  purposes.  If  the  situation  is  critical, 
that  is  near  category  limits,  then  other  means  of  analysis  must 
be  used. 

kef.  1  ''Instrument  Landing  Systems  Scattering"  Report 
No.  FAA-RD-72-137  (1972) 
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4.  TEST  CASE  FOR  THE  ILSLOC  COMPUTER  PROGRAM 


To  illustrate  how  the  computer  program  is  operated  a  very 
simple  test  case  (with  only  2  scatterers)  has  been  created  and 
run.  For  this  simulated  airport  the  program  computed  the  course 
width  as  4.01  degrees.  Both  anterna  arrays  were  set  at  an  eleva¬ 
tion  of  13  feet  above  the  ground  plane.  The  clearance  antenna 
array  was  used  as  the  origin  for  the  coordinate  system.  An 
80'xl00'x60'  hangar  and  75'xllO'  cylinder  were  placed  on  opposite 
sides  of  the  9,350  ft.  runway.  In  this  case  the  threshold  is 
10,000  ft.  from  the  course  antenna.  (See  illustration  -  Figure  2). 
Based  on  the  size  and  location.,  of  these  two  buildings,  the  model 
predicted  the  CDI  on  the  runwa)  centerline  and  for  a  clearance 
run  at  10,000  ft.  range. 

Using  this  model  for  input  values,  the  following  section 
presents  a  detailed  follow  through  of  the  main  program  steps. 

The  Mode  Card 

The  first  input  is  the  mode  card.  This  card  contains  informa¬ 
tion  on  the  type  of  localizer  antenna  used,  the  frequency  of  th 
ILS,  the  length  of  the  runway,  and  the  height  of  the  antenna. 
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The  card  format  is : 


Col.  Symbol 


Usage 


1-2  Mode 


*  1  (V-RING) 

=  2  (8-L00P) 

=  3  (WAVEGUIDE) 

=  4  (VACANT) 

=  5  (MEASURED  PATTERN)  indicates 

=  6  (MEASURED  CAPTURE  antenna 

EFFECT  PATTERNS)  type 

=  7  (THEORETICAL  PATTERN) 

=  8  (THEORETICAL  CAPTURE 
EFFECT  PATTERNS) 

=-l  (V-RING  CLEARANCE) 

--2  (8 -LOOP  CLEARANCE) 

=-3  (WAVEGUIDE  CLEARANCE) 

=-4  (MEASURED  CLEARANCE 
PATTERNS) 


11-20  FRQ 


Frequency  of  ILS  in  Mega  Hz 
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In  order  to  effectively  use  the  rest  of  the  mode  card 
columns  it  is  important  that  the  user  understand  the  coordinate 
system  used. 

The  x-axis  is  along  the  center  line  of  the  runway,  the 
threshold  being  in  the  positive  direction.  The  z-axis  is  vertical, 
positive  z  being  in  the  up  direction.  The  y-axis  completes  a 
right  handed  coordinate  system:  so  that  when  one  is  standing 
at  the  origin  facing  in  the  x-direction  positive  y  is  to  the 
left.  The  origin  is  used  as  a  reference  to  define  the  location 
of  scatterers,  antenna  system  components,  and  flight  path  sample 
points.  The  antennae  are  located  along  the  x-axis,  they  need  not 
be  at  the  origin;  as  in  our  test  case,  it  is  usually  convenient 
to  place  the  course  antenna  at  the  origin. 

Col.  Symbol  Usage 

21-30  XTH  Distance  from  the  origin  to 

the  threshold  of  the  runway, 
in  feet.  This  number  is  used 
for  both  flight  path  orientation 
and  for  course  width  determination. 
The  distance  is  given  in  feet. 

31-40  ZA (1 )  There  is  always  a  non-zero 

antenna  height,  and  it  is 
input  here. 

41-50  ZA(2)  This  will  be  the  clearance 

antenna  height  if  a  two  antenna 
system  is  used. 

Modes  1,  2,  and  3  provide  for  standard  localizer  antenna 
array  types.  These  antenna  arrays  are  predetermined,  the  only 
variable  being  course  width,  the  adjustment  of  which  is  controlled 
by  the  course  width  card. 

When  any  array  type  other  than  mode  1,  2,  or  3  is  used,  ad¬ 
ditional  antenna  array  description  cards  must  be  included.  Mode  5 
permits  the  input  of  a  measured  pattern  for  special  cases  on 
theoretical  studies.  When  this  mode  is  selected  additional 
pattern  cards  are  required.  One  pattern  card  must  be  used  for 
each  measurement.  The  angles  must  be  given  in  ascending  order. 

A  maximum  of  fifty  measurements  may  be  given;  if  less  than  fifty 
cards  are  used  a  termination  card  with  an  angle  greater  than  360 
degrees  must  be  inserted. 
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Format  of  Pattern  Card(s) 


Col. 

Symbol 

Usage 

1-10 

ANG 

Angle  of  measurement,  in  degrees 

11-20 

AFPP 

Amplitude  of  sideband  only  pattern, 
in  relative  units 

21-30 

AGPr 

Amplitude  of  carrier  plus  sideband 
pattern,  in  relative  units 

Mode  7  allows  the  generation  of  a  theoretical  array  pattern 
from  assumed  element  contributions.  The  antenna  is  to  be  a  linear 
array  of  elements  with  identical  radiation  patterns.  Each  element 
has  an  arbitrary  magnitude  and  phase  for  both  carrier  plus  sideband 
and  sideband  only  currents.  The  arrays  are  assumed  to  be  aligned 
parallel  to  the  y-axis.  All  elements  have  the  same  height,  as 
given  in  the  mode  card.  All  elements  have  the  same  x-coordinate  as 
given  on  the  course  width  card.  The  y-coordinate,  in  wavelengths, 
is  given  for  each  element  on  the  element  description  card.  There 
must  be  one  card  for  each  element  in  the  array,  to  a  maximum  of 
26  elements.  The  format  for  the  element  description  card  is: 


Col. 

Symbol 

Usage 

1-10 

DT 

Element  displacement  in  the 
y-direction  given  in  wavelengths 

11-20 

CT 

Carrier  plus  sideband  amplitude, 
in  relative  units 

21-30 

PC 

Carrier  plus  sideband  phase,  in 
degrees 

31-40 

ST 

Sideband  only  relative  amplitude 

41-50 

PS 

Sideband  only  phase,  in  degrees 

The  phase  of  the  sideband  only  currents  is  ideally  in  quadrature 
to  the  carrier  plus  the  sideband  currents.  This  9u  degree  shift 
is  added  by  the  program.  Thus  a  "PS"  inputted  as  zero  degrees  is 
internally  converted  to  90  degrees  out  of  phase  with  the  sideband 
portion  of  the  carrier  plus  sideband.  To  indicate  termination 
when  there  are  less  than  26  elements  used,  an  element  card  is 
placed  with  a  carrier  plus  sideband  phase  value  (PC)  of  more  than 
500. 


The  next  step  for  this  mode  must  be  the  input  of  the  horizontal 
radiation  pattern  for  the  individual  element.  This  pattern  will  be 
used  for  each  of  the  elements  previously  described.  The  input  is 
the  relative  signal  strength  measured  every  10°  starting  at  0  and 
proceeding  until  180°.  This  is  a  total  of  nineteen  amplitudes; 
the  values  are  read  in,  in  records  of  8F10.4  format,  for  a  total 
of  3  records.  This  gives  the  pattern  for  angles  from,  0°  to  180° 
and  since  the  pattern  is  assumed  to  be  symmetric  the  value  for  the 
negative  angle  will  be  the  same  s  a  positive  one  of  equal 
magnitude. 

There  are  two  methods  of  inputting  capture  effect  system 
descriptions.  The  most  general  way  is  to  input  each  antenna  array 
separately.  When  using  this  method  the  clearance  array  must  be 
input  first.  This  input  will  follow  the  same  steps  as  a  single 
array  system  except  that  the  mode  number  will  be  a  negative.  The 
negative  mode  card  and  the  pattern  or  element  cards  (if  any)  must 
be  followed  by  another  mode  card.  This  mode  for  the  course  array 
must  be  positive,  and  followed  by  the  necessary  pattern  or  element 
cards . 

There  are  two  cases  for  the  second  method  of  inputting  antenna 
array  descriptions.  The  first  case  is  used  if  both  course  and 
clearance  antenna  array  are  to  be  given  as  measured  patterns;  a 
single  mode  6  card  is  used  followed  by  two  sets  of  pattern  cards: 
the  first  set  is  for  the  course  antenna  array:  and  the  second 
set  for  the  clearance  antenna  array.  The  mode  6  is  converted  in¬ 
ternally  to  a  mode  5  for  each  array  and  these  values  will  appear 
in  the  output  listing.  In  the  second  case,  for  a  capture  effect 
system  which  uses  two  theoretical  arrays,  a  mode  8  is  used.  This 
card  is  followed  by  the  course  antenna  element  description  cards 
and  the  element  radiatiou  cards;  a  second  set  of  array  description 
cards  is  used  in  the  clearance  antenna.  As  in  the  mode  6  case, 
the  mode  8  is  converted  internally  to  two  mode  7's.  These  mode 
7's  will  appear  in  the  output  listing. 
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Figure  3.  Pattern  Card  Test  Case  Listing  (Cont.) 


The  label  card  follows  the  course  width  card.  This  card  is 
put  on  the  output  tape  ahead  of  the  CDI  records  for  this  flight. 
It  serves  as  an  identifying  record  and  is  the  label  placed  on  the 
graph.  Columns  1-80  are  used.  In  our  test  case  this  card  reads: 


I  THIS  IS  A  DEMONSTRATION  CASE  OF  STRAIGHT  LINE  FLIGHT. 

The  program  calculates  the  CDI  at  a  point  in  space:  for 
|  convenience,  the  program  will  permit  calculation  for  a  series 

|  of  points.  This  set  of  points  represents  samples  of  a  simulated 

|  flight  path. 

K 


|  The  program  allows  two  types  of  flight  paths.  A  straight 

I  line  flight  and  a  circular  orbit.  The  flight  path  card  has  one 

|  of  the  following  formats: 


Straight  Line  Flight 


1  Col. 

Symbol 

Usage 

1-10 

11-20 

l 

XMIN 

Starting  distance  from  origin, 
in  feet 

XMAX 

Ending  distance  from  origin,  in 
feet 

1  21-30 

a. 

DXR 

Spacing  between  sample  points, 
in  feet 

J  31-40 

g 

PHIR 

Angle  of  approach,  in  degrees 

1  41-50 

§ 

PSIR 

Glide  angle,  in  degrees 

|  61-70 

ZUP 

Height  of  aircraft  at  threshold 
in  feet 

XMIN  is  the  x-coordinate  of  the  starting  location  of  the 
aircraft  and  XMAX  is  the  x-coordinate  of  the  ending  location. 
The  sample  points  are  spaced  along  a  straight  line  so  that  the 
difference  in  x-coordinates  between  successive  samples  is  DXR. 
The  sign  of  the  DXR  will  be  set  by  the  program  so  that  the 
flight  goes  from  XMIN  to  XMAX  regardless  of  flight  direction. 

If  the  DXR  value  would  require  more  than  500  points  the  program 
will  adjust  the  magnitude  of  DXR  to  give  only  500  points.  In 
some  cases  a  flight  will  require  more  than  500  points.  If  this 
is  necessary  the  flight  must  be  broken  up  into  smaller  segments 


3j 


| 

3 

i 

3 


A 
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of  not  more  than  500  points  each.  The  procedure  for  doing  this 
is  explained  in  the  control  card  section.  The  flight  path  is 
oriented  in  space  so  that  an  extension  of  the  path  crosses  the 
threshold  at  the  altitude  of  ZUP  and  intersects  the  z-axis.  PHIR 
is  the  angle  between  the  flight  path  and  the  vertical  plane  through 
the  runway  centerline.  It  is  zero  for  a  flight  path  along  the 
centerline  of  the  runway  and  is  positive  for  an  incoming  flight 
(XMTN  greater  than  XMAX)  with  decreasing  y-displacement .  PSIR 
is  the  glide  angle  between  the  flight  path  and  the  horizontal 
plane.  It  is  zero  for  level  flight  and  positive  for  a  normal 
landing  approach.  The  flight  path  is  a  straight  line  as  de¬ 
scribed  above  except  when  the  x-component  is  less  than  XTH,  that 
is  if  the  aircraft  is  on  the  antenna  side  of  the  threshold.  In 
that  case  the  aircraft  altitude  will  be  set  up  to  ZUP. 

Thus  the  values  used  in  the  test  case  would  read: 


1-10 

40000. 

11-20 

20000. 

21-30 

-40. 

o 

1 

4 

0. 

41-50 

2.5 

51-60 

50. 

The  arc  flight  is  a  series  of  points  at  a  constant  height 
of  ZUP  and  at  a  constant  horizontal  distance  from  origin  of  R. 
MIND  is  the  starting  angle  for  the  arc,  that  is,  the  line  of 
sight  from  the  origin  to  the  point  makes  a  horizontal  angle  of 
MIND  degree  with  the  x-axis.  The  sample  points  are  spaced  at 
equal  angles  of  DXR  until  the  termination  angle  of  MIND  is 
reached.  As  in  the  straight  line  flight  the  sign  of  DXR  will 
be  adjusted  appropriately.  Likewise  the  magnitude  of  DXR  will 
be  set  to  yield  not  more  than  500  points.  Column  74  must  be 
set  to  1  to  indicate  a  circular  arc. 


| 


f 


1 


I 


■i 

A 


Circular  Orbit  Case 
Col.  Symbol  Usage 


I- 10  MIND 

II- 20  MAXD 


Starting  angle,  in  degrees 
Ending  angle,  in  degrees 
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Col. 

Symbol 

Usage 

21-30 

DXR 

Angular  spacing  between  samples 
in  degrees 

51-60 

R 

Radius  of  orbit,  in  feet 

61-70 

ZUP 

Height  of  orbit,  in  feet 

74 

ICF 

Must  be  set  to  1  to  indicate 
orbit  case 

Following  the  flight  path  card  must  be  the  velocity  card 
in  the  following  format: 


I  Col.  Symbol 

1-10  VEL 

I 


Usage 

Velocity  of  aircraft,  in  feet/sec. 
This  is  used  for  the  Doppler  Effect 
on  the  receiver.  The  sign  of  the 
velocity  will  be  made  to  agree  with 
the  directional  motion  from  DXR. 
Test  case  assumes  velocity  of 
200  ft. /sec. 


At  this  point  we  have  described  the  antenna  system  and  the 
trajectory  of  the  aircraft;  the  derogating  surfaces  in  proximity 
to  the  ILS  must  now  be  described.  The  program  will  simulate 
scattering  from  rectangular  or  cylindrical  surfaces.  We  will  now 
describe  the  method  of  inputting  scatterers  to  simulate  derogating 


structures . 


The  next  card  describes  either  the  scatterer(s)  or  output  and 
control.  The  usage  is  determined  by  the  value  of  the  ID  field 
in  columns  1  to  2.  An  ID  of  -1,  1  or  2  is  used  for  scatterers, 
while  the  other  values  are  used  for  control. 


i 

% 

Col. 

Symbol 

Usage 

■-4 

1S-20 

XW(3; 

Z-coordinate 

; 

3 

26-30 

ALPHA 

Angle  between  base  and  x-axis,  in 

i 

degrees 

- 

* 

31-35 

DELTA 

Angle  of  tilt,  in  degrees 

= 

36-45 

WW 

Width  of  rectangle,  in  feet 

i 

5 

46-55 

HW 

Height  along  rectangle,  in  feet 

g 

The  scatterer  is  a  rectangle  with  the  reference  point  at  the 
middle  of  the  base.  The  rectangle  is  assumed  to  be  of  infinite 
conductivity  and  zero  thickness.  It  also  has  only  one  side.  This 
can  be  thought  of  as  the  front  surface  of  a  metal  wall.  A  wall 
with  zero  x-,  y-,  and  z  coordinates  and  an  alpha  of  zero  is  located 
at  the  origin  with  surface  of  the  wall  facing  in  the  negative  y 
direction  (Figure  4,  case  I).  A  positive  increase  in  alpha  rotates 
the  wall  about  the  z-axis  in  a  counterclockwise  direction  when 
viewed  from  above.  Thus  an  alpha  of  ninety  degrees  faces  the  wall 
in  the  positive  x  direction  (Figure  4,  case  II).  Alpha  is  the 
angle  between  the  vertical  projection  of  the  base  of  the  wall  in 
the  xy-plane  and  the  x-axis,  measured  in  degrees.  Delta  is  the 
angle  between  the  surface  of  the  wall  and  the  vertical  direction, 
in  degrees.  A  delta  of  zero  is  a  wall  perpendicular  to  the  ground 
and  a  decrease  in  delta  rotates  the  wall  about  the  baseline  in  a 
direction  so  that  a  delta  of  minus  ninety  is  a  horizontal  wall 
facing  down  (Figure  4,  case  III).  WW  is  the  width,  in  feet,  of 
the  wall  measured  along  its  base  and  HW  is  the  height  measured 
along  the  surface  at  right  angles  to  the  base.  If  the  wall  is 
oriented  in  such  a  fashion  that  the  line  of  sight  from  the  antenna 
to  the  wall  passes  through  the  back  and  not  the  front  of  the  wall, 
the  program  will  ignore  the  wall  in  the  simulation. 

An  ID  of  *1  is  used  with  the  above  format  to  describe  a 
negative  wall.  This  ID  is  used,  for  example,  to  create  a  wall 
with  a  rectangular  hole  in  it.  The  entire  surface  is  used;  the 
hole  is  then  subtracted  by  inputting  a  second  card  with  an  ID  of 
-1  and  the  size,  location,  and  orientation  of  the  hole. 


An 

ID  of  2  is  used  for  a 

cylindrical  scatterer  with  the 

t 

following  format: 

-= 

Col. 

Symbol 

Usage 

- 

1-2 

ID 

Must 

be  a  2 

1 

3-8 

XW(1) 

1 

1 

9-14 

XW(2) 

y-  1 
| 

i  coordinates  of  the 

r  reference  point,  in  feet 

4 

15-20 

XW(3) 

2- 

1 

3 

f 

m 

36-45 

HW 

Diameter  of  cylinder,  in  feet 

i 

46-55 

HW 

Height  of  cylinder,  in  feet 

5 

The  reference  point  is  located  at  the  base  of  the  cylinder  on 
the  axis  of  rotation  of  the  cylinder.  The  diameter  is  W*f  feet,  with 
the  base  parallel  to  the  xy  plane  at  an  altitude  of  XW(3)  feet.  The 
cylinder  extends  upward  for  HW  feet  with  the  axis  of  rotation  in  the 
vertical  direction.  The  cylinder  is  assumed  to  have  infinite 
conductivity. 

After  an  ID  of  -1,  1  or  2,  the  program  will  calculate  the 
electric  field  at  the  surface  of  the  scatterer.  This  will  be 
calculated  from  the  signal  from  the  transmission  antenna  array 
and  from  the  ground  reflection  of  the  transmitted  signal.  Then, 
for  each  receiver  point  along  the  flight  path,  the  program  will 
calculate  the  electric  field  at  that  location  from  the  scattered 
signal:  from  both  the  scatterer  and  reflected  from  the  ground. 

Thus,  the  signal  is  received  from  four  paths:  transmission 
antenna  to  scatterer  to  receiver;  antenna  to  ground  to  scatterer 
to  receiver;  antenna  to  scatterer  to  ground  to  receiver;  and 
antenna  to  ground  to  scatterer  to  ground  to  receiver.  This 
signal  is  decomposed  into  co^lex  components  induced  in  the 
receiving  antenna  at  the  different  carrier  and  sideband  fre¬ 
quencies.  The  program  then  .’oops  back  to  read  in  another  ID 
card,  permitting  the  summation  of  the  effects  of  many  scanterers. 
This  allows  the  simulation  of  complex  structures  by  breaking 
them  up  into  cylinders  and  rectangles. 
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In  the  test  case,  we  have  only  inputted  three  scattering 


surfaces.  This  was 

and  the  cylinder  are 

cards  read: 

done  because  only  two  sides  of  the  hangar 
illuminated.  The  values  for  the  scatterer 

Col. 

First  card 

Second  card 

Third  card 

1-2 

1 

1 

2 

3-8 

6000. 

5950. 

7500. 

9-14 

1100. 

1130. 

-1000. 

15-20 

0. 

0. 

0. 

26-30 

10. 

-80. 

0 

31-35 

36-45 

100. 

60. 

75. 

46-55 

80. 

80. 

110. 

After  all  the  scattered  have  been  input,  a  control  card  is 
inserted  to  terminate  the  run.  The  control  card  format  is: 

Col. 

Symbol 

Usage 

1-2 

ID 

not  -1,  1, 

or  2 

When 

a  control  card 

is  read  in,  the  program 

will  add  the  direct, 

a>:d  ground  reflected  signal  from  the  transmission  antenna  to  the 
scattered  signal  summations,  thus  giving  the  total  received  signal. 
The  program  then  calculates  the  CDI  that  would  be  seen  at  each  re¬ 
ceiver  point,  and  outputs  the  label,  a  ht - ier  record  describing 
the  flight  path  and  the  values  of  the  CDI  on  output  tape.  If  the 
ID  is  equal  to  zero  the  program  also  outputs  additional  records  for 
the  strengths  of  sideband  and  carrier  signals  from  course  and 
clearance  (if  any)  antenna  arrays.  The  field  summations  are  then 
cleared  for  the  next  run. 

The  program,  having  finished  the  previous  run,  now  proceeds 
with  the  next  input.  The  next  run  is  generated  by  looping  back 
to  a  point  in  the  input  stream,  determined  by  the  value  on  the 
control  card. 
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Once  an  input  sequence  has  begun  the  inputs  following  in  the 
standard  order  must  be  given.  The  user  must  also  keep  in  mind  that 
all  values  on  cards  given  before  that  entry  point,  in  the  previous 
run  are  still  in  effect.  The  standard  order  is: 

MODE  CARD 

(measured  pattern  for  modes  5  and  6  or  current 
description  for  modes  7  and  8) 

(second  mode  card  and. patterns  of  currents  if 
first  mode  was  negative) 

COURSE  WIDTH  CARD 

LABEL  CARD 

FLIGHT  PATH  CARD 

VELOCITY  CARD 

(set  of  scatterer  cards) 

CONTROL  CARD 

The  value  of  the  ID  on  the  control  card  guides  the  looping 
in  the  following  manner: 

Value  of  ID  Next  card  to  be  read  in 

0  MODE 


3-10 

SCATTERER 

11-15 

LABEL 

16-20 

MODE 

21-50 

COURSE  WIDTH 

>50 

WILL  CAUSE  THE  PROGRAM  TO 

TERMINATE  AFTER  OUTPUTTING 
THE  LAST  CDI 


The  looping  permits  the  repetition  of  a  run  with  changes  in 
some  or  all  of  the  variables.  For  example,  ID  values  3  through  10 
permit  a  run  with  the  same  antenna  system  and  flight  path  as  the 
previous  case,  but  with  a  new  set  of  scatterer  inputs. 

ID  values  11  -  15  permit  a  new  flight  path  description  and 
scatterer  set  to  be  input.  This  looping  method  can  also  be  used 
for  flights  that  would  require  more  than  500  points.  For  reliable 
simulation,  the  spacing  between  receiver  points  (DXR)  should  be 
small  enough  so  that  the  change  in  CDI  between  successive  points 
is  not  more  than  ~ 20%  of  the  peak  value.  Thus  for  long  flights  the 
flight  path  must  be  broken  up  into  shorter  segments.  If  the  number 
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of  segments  of  this  path  does  not  exceed  4,  the  plotting  program  will 
connect  them  on  a  single  graph.  The  control  for  this  joining  is  the 
ID  number.  If  the  flight  path  finishes  with  an  ID  of  11  -  13,  the 
graph  of  the  next  flight  will  continue  the  line  of  the  graph.  A 
long  flight  may  be  broken  up  into  as  many  as  four  segments:  with 
three  segments  terminating  in  11  -  13  and  a  fourth,  and  final  seg¬ 
ment,  terminating  in  14  or  15.  The  flight  segments  must  appear  in 
the  order  in  which  they  are  to  be  flown,  so  that  the  XMIN  of  one 
section  is  the  XMAX  of  the  previous  section.  For  each  segment 
the  programmer  must  re-input  the  same  scatterors.  If  only  one 
segment  is  to  be  plotted  the  control  card  should  read  14  or  15. 

ID's  16  through  20  start  inputting  at  the  mode  card,  thus 
allowing  a  completely  new  run. 

An  ID  of  21  through  50  uses  the  same  antenna  description,  but 
starts  the  inputting  at  the  course  width  card.  This  permits  the 
course  width,  clearance  strength  and  antenna  location  to  be 
varied. 

The  program  is  terminated  after  an  ID  greater  than  50  is  en¬ 
countered.  The  direct  signal  will  be  added,  and  the  CDI  will  be 
outputted  before  the  program  stops.  The  program  will  also  stop 
if  an  end-of-file  is  encountered  while  the  program  is  attempting 
to  read  any  input  card,  or  if  certain  of  the  variables  are  of  im¬ 
proper  value.  In  these  cases  the  program  terminates  immediately, 
without  outputting  the  last  case. 

The  input  of  the  test  case  flight  path  was  done  in  four 
segments.  The  first  segment  is  from  40,000'  to  20,000',  the 
second  segment  is  from  20,000'  to  12,500',  the  third  segment  is 
from  12,500'  to  11,000'  and  the  last  is  from  11,000*  to  10,000'. 

An  additional  case  for  a  simulated  clearance  flight  by  a  circular 
orbit  has  also  been  included.  The  input  cards  for  these  test 
case  flights  are  shown  in  Figure  5. 
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THIS  Is  A  DEMONSTRATION  CASE  OF  STRAIGHT  LlNt  FLIGHT 
Aooon.  ?noon,  -An.  , 

>nn.  '■  3 

isnnn,  lino.  10.  100.  80. 

1],n*  -80.  60.  gO. 

27500.  -1000.  0.  0.  0.  75.  110. 

t? 

THIS  IS  A  DEMONSTRATION  CASE  OF  STRAIGHT  LlNt  FLIGHT 
20ono.  12500.  -15.  p.s 

>00. 

16000.  1100.  10.  100.  80. 

1  soSn.  11V.  -80 .  60.  80. 

27SOrt.  —loon.  o.  o.  0.  75.  110. 

s 

THIS  IS  A  DEMONSTRATION  CASE  OF  STRAIGHT  LINE  FLIGHT 
12500.  1  loon.  -3. 

0. 

16ooo.  lino.  10.  100.  80. 


I 


16950.  1130. 
27500.  -1000a  O. 
5 


60. 

0.  75. 


80. 

110. 


THIS  IS  A  DEMONSTRATION  CASE  OF  STRAIGHT  LINE  FLIGHT 

11000.  10000.  -2.  2.5 

no.  * 


lfioon.  lino, 
16050.  1130. 
2760o.  -loon.  o. 
6 


100. 

60. 

0.  75. 


80. 

SO. 

110. 


THIS  IS  ORRIT  CASE  WITH  SIGNAL  STENGTHS 
180.  180.  0.72 

200. 

16000.  1100.  10.  100.  60. 

15950.  1130.  -80.  60.  80. 

27500.  -1000.  0.  0.  0.  75.  110« 


10000, 


Figure  S.  Flight  Case  Inputs 
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EFN  SOURCE  STATEMENT 
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US  SIK’SLE  REFLECTION  INTERFERENCE  PROGRAM  ILSLOC 

this  program  simulates  the  effects  of  rectangular 

AND  CYLINDRICAL  SCATTERERS  ON  THE  lOCAL!*S»  PART 

OF  THE  ILS.  THESE  COMMENTS  SERVE  AS  THE  PROGRAM  DESCRIPTION 

FOR  THE  USER.  A  USER.S  MANUAL  HAS  SEEN  WRITTEN  AND 

THIS  COMMENTARY  IS  WRITTEN  ASSUMING  THE  USER  HAS  READ  IT. 


ILBL  IS  USED  TO  IDENTIFY  THE  SIGNAL  STREAM  OUTPUTS  AS 
TO  TYPE  AND  SOURCE.  THE  FIRST  CHARACTER  IS  ,S#  FOR 
SIDEBAND  ONLY  SIGNALS  OR  ,C,  FOR  CARRIER  PLUS  SIDEBAND. 

the  second  pair  are  .cr,  for  course  antenna  op  »cl.  for 

CLEARANCE. 

DIMENSION  ILBL(S) 

data  ILBL/4HC  CR.4HS  CR.4HC  CL.4HS  CL.4H  001/ 


LOGICAL  EOF 

OIMENSION  MEMO(14),DF(501) 

COMPLEX  9ESF.FAC.CE 

COMPLEX  £P,eE.EM,EC,2EM),?0<i>.EWR,rPP.GPPiFPM,GPM» 

2  CS(29,2) ,S0(25,2) 

COMPLEX  |JM,2JP,ZJPC(2)»2JMC(2> 

COMPLEX  ZP(900).ZPe(502.2>.2M(5!?l0),?Me(500.2> 

DIMENSION  XXRY(500,4> 

DIMENSION  VCD (900. 2) .VPO (900. 2) . VMM (50?, 2) 

DIMENSION  XW(3).Xtf0(3> 

DIMENSION  AN(3) 

DIMENSION  AF00(9).PHS(9) 

DIMENSION  XY(13) 

REAL  LAMBDA 

COMMON/CD/  ARaO(50).AFPP(50).ASPP(901.BRAO(90>.BPPP(00)»BGPP(901 

COMMON  /AB/  2JM.2JP.ZJPC.1JMC 

COMMON  2P.ZPC.2M.ZMC. VCD. VPO. VMD 

COMMON  /VAR/  SM.SNCUT,SNCUD.SNCUC(2)#VPC(2)»VMC(2I 

common  /SU9/  MODE, ICP, FRO. LAMBDA, pi, RAOO. PHI <S).P5I(S),NEL.XTh, 

1  XXA(3),YA,ZA(3),RA(3) 

COMMON  /ANT/  LOC.FPP.FPM, GPP.GPM, EWR(4, 4) , CWA (2) » AS»CLS»0E(29,2i , 
.  CS,SO.ET(20.2>.NO(2) 

EQUIVALENCE  (ZP(1), 20(1)), (XXPY(1, 1)i ZP(&)) 

OATA  RAO/57. 2997799/ 

CP  AND-  CM  ARE  THE  AMOUNTS  OF  MODULATION  ON  THE  CARRIER 

for  the  carrier  plus  sideband,  cp  is  the  course  modulation 

AND  CM  THE  CLEARANCE. 

DATA  CP, CM/. 2,. 2/ 


THE  OUPUT  of  the  SIMULATION  is  ON  unit  8,  A  TAPE  WITH 
WRITE  RING  SHOULD  BE  PLACED  THEREON. 
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£Fm  SOURCE  STATEMENT 
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PrM'-.r  9 


Nr  is  the  count  of  the  case  being  simulates  tt.s  value  ic  written 

ON  TWF  TAPE  WITH  THE  OUTPUT  RECORD.  THIS  WOULD  ALLOW 
SEARC'-HG  POP  A  PARTICULAR  case  BY  NUMBE*. 

'Ksg! 


this  IS  Tu£  STARTING  POINT  FOR  A  SIMULATION.  IT  je  ALSO 
ENTERED  Frs  A  PESTaRT  FOLLOWING  an  ip  of  e  0*  16  TO  20. 

1  CONTINUE 

-ode  * 

NEL  IS  THr  NUMBER  OF  ANTENNAE  IN  TmE  SYSTEM.  DEFAULT 
CCNOI TI0N  IS  ONE  ANTENNA 

N£L  *  1 


EWR  IS  A  OOMP'tX  MATRIX  CONTAINING  THE  SIOEBAND  ELECTRIC  FIELD 
DESCRIPTION  PRODUCED  BY  Tuf  ANTENNA  SUBROUTINE.  EwR(I.J) 

IS  THE  FIELD  FOR  THE  ,I,Th  ANTENNA.  AND  THE  »J»  VALUES 
HAVE  THE  pOLLOWING  SIGNIFICANCE! 

J  USAGE 

1  SIDEBAND  PORTION  OF  CAPRIER  PLUS  SIOEBAND 

FOR  THE  COURSE  SECTION  OF  THIS  ANTENNA 

2  SIOEBAND  ONLY  FOR  THE  COURSE 

3  SIOEBAND  PORTION  OF  CARRIER  PLUS  SIOEBAND 

FOR  THE  CLEARANCE  SECTION 

4  SIOEBAND  ONLY  FOR  THE  CLEARANCE 

THIS  SUBROUTINE  CALL  IS  USED  TO  CLEAR  EWR  BEFORE 
STARTING  THE  SIMULATION 

CALL  CLEAR(EWR,16) 


THIS  IS  A  TEST  FOR  ENO-OF-FILE  ON  CARD  INPUT.  THE  CALL  TO 
EOF  ARMS  the  JNTER'JPT,  AT  END  OF  FILE  ON  UNIT  5  I MTrRUPT  IS 
TO  STATEMENT  55. 

2  CONTINUE 

IFCE0F<5>>  GO  TO  58 


THIS  IS  THE  INPUT  POR  THE  MODE  CARP.  THE  VAR! ARLES  HAVE 
THE  FOLLOWING  USESt 

SYMBOL  USE 

mode  antenna  type 

*1  V-RING  COURSE 

*2  8-LOOP  COURSE 
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HAJM  -  EF.M  SOURCE  RTATfMCNT  -  IPNCS5  - 

*3  WAVEGUIDE  course 

■4  NOT  USED 

•3  measureo  course  pattern 

■6  MEASURED  COURSE  AND  CLEARANCE  PATTERNS 

•7  theoretcal  course  array 

■3  THEORETICAL  course  *no  clearance  array 

••1  V-RING  CLEARANCE 

•-?  8-LOOP  CLEARANCE 

•-1  WAVEGUIDE  clearance 

•-3  MEASURED  CLEARANCE  PATTERN 

■-7  THEORETICAL  CLEARANCE  ARRAY 

ERQ  FREQUENCY  OF  TRANSMISSION 

XTH  DISTANCE  to  threshold 

2A  ( I )  .  I# TH  ANTENNA  HEIGHT 

ORIGIN  IS  AT  THE  CENTER  OF  COORDINATE  SYSTEM. 

X-AXIS  IS  ALONC  RUNWAY 

2-AXIS  IS  STRAIGHT  UP 

Y-AXIS  COMPLETES  A  RIGHT  HANDED  SYSTEM 

RFAO  (5,lfl01>  MOOE.FRQ.XTH.2A 

this  is  a  test  for  invalid  antenna  type,  the  program  aborts  in  case 

OF  ERROR.  THIS  IS  USUALLY  CAUSED  RY  OMISSION  OF  OTHER  CARDS 

which  cause  something  other  than  a  mode  card  to  be  Reao  at 

THIS  POINT. 

IF  t  MODE  .GT.  8  )  GO  TO  38 
IF«  MODE  .LT.  -7)  GO  TO  38 
ir(  MSCE  ,E3.  0)  GO  To  98 

this  is  test  for  negative  mode  indicating  clearance  antenna, 
if  mode  is  positive  flow  is  to  statement  4 

IF C  MODE  .GT.  0  )  GO  TO  4 

ICP  is  the  .antenna  type  for  the  clearance  antenna 
leP  •  -  MODE 

if  there  is  a  clearance  antenna  then  the  numbe®  of  antennae 

IS  SET  TO  2 . 

NEL  •  2 

IF  the  CLe-ARENCE  antenna  IS  SPECIFIED  BY  A  MEASURED  PATTERN  IT  IS 
NOW  READ  in  by  subroutine  pattrn. 

IFC  ICP  .EQ.  3  )  CALL  PaTTRN(BRAD.PFPP.BGPP) 


rtiwpiw^ 
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C  I '  TMC  CLfASC.vce  Awtr.MMA  IS  SPECIFIED  a?  A»RAY  PARAMETERS  THE  INPUT 
C  T'Ta  FOP  TH"  iRSAY  IS  NOy  READ  J\*  py  C?RNT?. 

c 

IFU'V  .EO.  7)  CALL  CRRNTS  (DC ( 1  * 2) » 05 ( i» 2)  * SQ<  1 , 2 )  f  ET(1, 2 > . NO (2) ) 
C 
c 

c  Tu;  ~!.0U  is  <0W  SACK  TO  STATEMENT  7  TO  READ  IN 
C  .•r.VE  ?4S0  F OP  COURSE  ANTEVA. 

C 

S'*  TO  ?. 

C 

C 

C  TMIS  *S  The  INPUT  SECTION  FOR  THE  COURSE  ANTENNA  IF  PATTERNS  OR 
C  ARRAY  OfRrDI »TiON  MUST  9S  5 1  YEN*  OTHERWISE  FLOW  IS  To  THE 

c  initialisation  section. 

c 

4  ire  m:.ce  .lT.  s  j  go  to  a 
c 
c 

c  this  staTehent  controls  the  input  method,  pattern  or  array, 

C  ACCORDING  to  mode  TYPE, 

c 

IF  (MODE  ,rT,  6)  GO  to  5 
CALL  PiTT3;s{ARAD,AF0P.  AGPP) 

C 

C 

c  this  is  t?  input  t«e  second  pattern  fop  clearance  antenna  if 
C  ''ODE  ts  6. 
c 

IF (  MODE  .EQ.  5)  GO  TO  6 
CALL  PATTRN(PRAD.BFPP,9GPP> 

C 

C  Thf  NUMBER  of  ANTENNAE  and  the  ICP  type  ARE  SET.  TuE»>  FLOW  IS  TC 

c  initilaiiation. 
c 

NEL  »  2 
MODE"? 
irP  •  5 
GO  TO  6 
C 
C 

c  this  is  the  input  eor  coupse  array  DATA. 

S- CALL  'crrnts'  foEVcsVso',  FT, NO (1  > } 
c 
c 

C  this  TEST  IS  FOR  CLEARANCE  ARRAY  IF  MODE 
C  IS  TYPE  B 
C 

IF  (  MODE  .EO.  7)  GO  TO  6 

CALL  CPRNTS  <D£(1.2).CS(1.2).S0(1.2)iET(1,2).ND<?)) 

M0DE«7 

1CP»7 

NEL«2 

C 

C 
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efn  scupc*  STATEMENT 


IFh(S> 


C  THIS  IS  ThF  IMTULI?ATinM  SECTION.  LAM°DA  IS  TH£  WAVELENGTH 
C  IK  FEET  4-D  A  *  IS  THE  PHASE  SHIFT/.-  IRTaNOE  fN  »ADIa*JS/F00T. 

C  YA  IS  ThE  Y-COOPfllNATE  GF  THE  ANTEvNAE.  Tu»S  IS  ASSUMED  TO 
C  PE  IE°0  I «  ALL  CASES. 

C 

6  LAMpn4sii<!?0./Fqa/12. 

AK*2.*?I/L4MP0A 

YA*C!.n 

*2  *  r. 

IF(E0**(5>)  GO  TO  5« 


THIS  IS  The  course  WIDTH  INPUT. 

XXA(l)  is  THE  X-COCRDR'ATE  OF  THE  COURSE  AmT£n\a 
XXA<2>  IS  T HE  X-COOROINATE  OF  THE  CLEARANCE  ANTENNA 
cm  is  the  course  width 

CLS  IS  THE  ratio  CF  CLEARANCE  TO  COURSE  SIGNAL  strength, 
READ  <*,iP0?>  VXA.CW.CLS 


set  THE  PEFaJlT  CONDITION  ON  CLS  OF  i. 
Irc  CLS  .LE.  f!.a  )  CLS  •  1.3 


CWACI)  IS  THE  COURSE  width  ADJUSTMENT  ON  The  ,i,Th  antenna 
IT  SETS  The  Sideband  TO  CARRIER  RATIO.  THE  CLEARANCE  ANTENNA 

ccwa(2>)  is  always  1.3  .  the  course  width  is  adjusteo 
8v  varing  the  course  antenna  <ew4(t>>. 

CWA(l)  a  l.e 
CWA(2)*1»0 


THE  PROGRAM  will  NOW  CALCULATF  the  COI  FOR  2.5  DEGREE  COURSE 

OFFSET,  this  is  USED  TO  NORMALIZE  THE  SIDEPAjC  level  to 

ACHIEVE  The  desired  course  WIDTH.  LOC  IS  the  type  OF  antenna 

used  by  the  antenna  subroutine,  psicd  is  the  angular  altitude 

of  THE  REFERENCE  POINT  and  PHKl)  IS  ThE  A2l“UTH  OF  THE  POINT. 

Psm>  «  5.E-73 
PH  1(1)  a  ?.5»PAD3 
LOC  •  MODE 

c 

c 

C  THE  MODE  IS  USED  TO  DETERMINE  WHICH  ANTENNA  SUBROUTINE  TO  CALL. 

c  csp  is  the  standard  antenna  routine,  it  covers  the  v-ring 

C  8-LOOP  AND.  WAVEGUIDE.  LNaR  IS  THE  ARRAY  ANTENNA  SUBROUTINE. 

C  ANTP  IS  ThE  MEASURED  PATTERN  SUBROUTINE.  THE  SUBROUTINE  WILL 
C  RETURN  FP°  ANO  GPP  FOR  Tw£  POINT  AT  PHI,®$!  AND  UNIT  RANGE. 

C  FPP  IS  THE  SIDEBAND  ONLY  LEVEL.  GPP  IS  THE  SIDEBAND  LEVEL 
C  FOR  The  Carrier  PLUS  SIDEBAND.  AFTER  THE  peTURN,  flow  IS  TO 
C  STATEMENT  9. 

C 

IF  (MnrE  .GE.  7}  GO  To  8 
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«tiv  -  £fw  source  statement  -  irN(S)  - 

IFfMOpE  .GE.  5)  GO  TO  7 
CALL  CSP 
50  TO  9 

7  CALL  aVT»  (FPP,GpP,ARAO.AFPP,AGPP) 

50  TO  9 

ft  CALL  !>-AP  (  FPP.GPP.PHl.OE.CS.SC.ET.ND) 

GO  TO  9 
C 
C 

C  THE  StGNAl  LEVELS  ARE  IN  FPP  AND  GPP.  TEMP  !S  THE  APPARENT 
C  COURSE  WIDTH  WITH  CHA.S  OP  i.E. 

C 

9  T£HP«  1. 9375/REAL (PPP/GPP) 

C 

C 

C  THE  COURSE  WIDTH  READ  in  is  USEO  IF  IT  JS  LARGE*  than  3  DEGREES 

c  otherwise  the  standard  value  bt  faa  specifications  is 

c  DETERMINED  ANC  THIS  VALUE  USED,  Tw£  COU*SE  MIRTH  IS  LIMITED 
C  TO  A  RANG?  or  3  to  6  DEGREES. 

C 

if (  cw  -  s. a  )  10.10.11 

12  CW  ■  2,*ATaN(350./XTH  )  •  RAO 
IF(  CW  ,LT.  3.3  )  CW  B  3.? 

IFICW  .GT.  6c0)  CWbS.0 
C 
C 

C  THE  CWA(l)  IS  ADJUSTED  TO  PRODUCE  THE  DESIRED  COURSE  WIDTH. 

C 

11  CWA(l)  s  temp/cw 
C 
C 

c  the  values,  reaq  in  and  calculated,  for  the  antenna  sy$tem<s> 

C  ARE  OUTPUT  on  THE  LINE  PRINTER  (ASSUMED  TO  BC  UNIT  6) 

C 

WRITE (3. If 03)  MODE. ICP, FRO, XTW.2A.XXA.CW 
WR1TE(6,1000 )  TEMP.CWA 
WRITE (A. 1002)  CLS 
C 
C 

c  this  is  The  loop  back  point  for  new  FLIGHT  PATH,  ID, S  11  TO  15. 

C  MEMO  IS  The  label. FOR  HEADER  RECORDS  and  GRAPHS. 

C  INPUT  DATA  FOR  FLIGHT  PATh» 


C 

XM'IN 

STARTING  POINT 

c 

xmax 

ENOING  POINT 

c 

oxp 

SAMPLE  POINT  SPACING 

c 

PHIR 

ANGLE  OF  APPROACH 

c 

PSIR 

GLIDE  ANGLE 

c 

R 

RAOIUS  OF  ORBIT 

c 

ZUP 

ALTITUDE  AT  THRESHOLD  OR  OF  ORBIT 

c 

ICF 

FLAG  0  FOR  STRAIGHT  LINE,  1  FOR  Oi 

C 

14  CONTINUE 

READ  (5,1105)  MEMO 
WRITEC6.10R4)  MEMO 

READ  (5, 1IS4)  XMIN. XHaX, DXR.PHIR.PS I R.R.IUP. ICF 
C 
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c 

c 

c 

c 

c 


the  sion  -f  oxr  is  adjusted  for  plight  from  shin  to  xmax. 

OVR«SlCN(*:yR#  (XHAX-XMrj) } 


2  THP  VPLOCITY  OF  THE  AIRCRAFT  IS  INPUT. 


c 

c 

c 

e 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 


c 

e 

c 

e 

•  C 


R£A0  (S,i;?6)  vel 

WRITE  (6.1727)  VEL 

the  sign  of  the  velocity  is  set  to  agbee  with  that  of  dxr. 
VFL«S:GN<"EL.3XR) 

The  NliMaEo  oe  RECEIVER  POp'TS  IS  DETERMINED.  IF  ThI*  JS 
LESS  than  *02  FLOW  PROCEEDS  TO  STATEMENT  14.  OTHERWISE  THE 
MACNITOOE  OF  DXR  IS  INCREASED  TO  GIVE  DULY  5*1  POINTS, 

NR  slFlXC  ( XHAX-XW IM ) / jXR  *  1.) 

IFCNR  .LT.  1>  GO  TO  15 
lFtNR-501)  14,10.15 

15  WDITE  (0,1?»P) 

DXR«(vMAX-XHIN)/502. 

IF(aSoCOXR)  .LT.  l.E-5)  GO  TO  «S 
NP  ■  5»l 

16  CONTINUE 

the  flight  path  description  is  output,  t»s  format  3fing  detcrmincd 

by  THE  TYPE  OF  PLIGHT.  IN  THE  CASE  OF  STRAIGHT  LINE  THE 
NECESSARY  CONSTANTS  pqR  DOPPLER  EFFECTS  AND  POSITION  ARE 

determined, 

AFTER  OUTPUT  FLOW' IS  TO  STATEMENT  19. 

IF  CICF)  15,14,17 

17  WRITE  (6.1015)  XMIN,XmaX,OXR,XTH,ZUP. ICF 
GO  TO  19 

18  CONTINUE 

WRITE  '(6.1019)  XMIN,Xmax,3XR,Phir,p$jr,xTw,ZUP 

phimphir/rad 

P$tMPS!»/RAO 

SPSI  >  SIN(PSIR) 

TaNS**SPSI/COS(PSIR) 

TaN#R*SIN(PHIR)/COS(PhIR) 

V.X«VEl«COS(PSIR)»COS(PhIR) 

VY«VEl«COS(PSIR)*EIN(pui») 

V?»VEl*SIN(PSIR) 

19  CONTINUE 

these  constants  are  filter  factors  for  the  assumed  modulation 

FILTERS, 

Ta«1./15. 
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efn  source  statement 


ITn(S) 


t4/: 


»?4  S  PI»Ta 
Wf*T  a  6’.*«»TA 

We>'T  =  r:.#pta 

W1*\?T  =  1=5?,«PTA 


THIS  *  **  T-.it  L ’OP  BACK  POINT  TO  START  A  NEW  SIMULATION  WITH 
0S*EVI‘.i»3  4*T£'*NA  SYSTEM  a  ns  FlIGHT  PATH,  the  complex  fielo 
SUKMATiO"  MATMCI£S  ape  CLEAPEP,  The  Case  NUMRrR  is 
INCREMENTED  V  ONE  AMO  THE  LlN'EPRMCR  HEADERS  ARE  WRITTEN. 

??  CONTINUE 

CALL  rLt^(lP.43<H2) 
mc  «  *jc  ♦  i 
NfclTE  <6*1«P) 

WRITE  (6.1211) 


this  is  The  input  for  a  new  scatteber  or  control  card,  the 

PT.RMAT  ANn  USAGE  OF  THE  VARIABLES  WILL  BE  FOUNC  IN  TmE  USER.S  MANUAL. 
21  READ  (3,1"12>  10. XW( 1 ) , XW(2) »XW (3) « ALPHA. DELTA* WW.HW 


A  NEGaTi Vr  i:  is  USED  ON  A  SCATTEBER  TO  CAUSE  THE  FlELOS  TO 
BE  SUBTRACTED  PROM  THE  SUM.  THUS  |OA  IS  USCO  TO  DETERMINE 

the  type  op  s~atterer  and  the  sign  of  io  is  used  f?r  the 
sign  ''Termination  of  the  fields. 

!QAalAnS( TO) 


The  receiver  POINT  LOCATION  VARIABLES  arc  initialized,  xr  IS 
The  X-CGORDINATE  of  the  location,  zr  is  the  i-cooroinate 
AND  C0E6  ?s  the  azimuth,  the  use  of  these  VARIABLES  IS  CONTROLLED 
BY  THE  VALUE  OF  ICF. 

IF  (ICF)  23.23.Z2 

22  C0EG»XM!N-OXR 
ZRRfUP 

GO  TO  24 

23  XP*XKIN-OXR 

24  CONTINUE 


IF  IDA  is  not  i  OR  2  THEN  THIS  CARO  IS  A  CONTROL  CaRQ  AND 
FLOW  PASSES  TG  STATEMENT  43  TO  OUTPUT  THC  CO!  AND  FOR 
LOOPING  CONTROL* 

IFUDa  .CT.  2)  60  TO  43 
IFUOa  ,E3.  B)  60  TO  43 


xv  is  an  array  of  data  on, the. antenna  ano  FlICHT  Path  and  If 
output  as  part  of  the  header  rccoro  on  the  output  tape. 


34 


-  £Fn  SOURCE  STATEMENT  -  IFNIS)  - 


Kv(jj)  a  a  ,i« 

XY(2>  •  »5lR 
XYCJ)  *  2'JP 
XYU)  *  pLOAT(MC) 
XV(5J  a  V£L 
YYt#>  s  P'.CAT (MODE ) 
XY(7)  *  "U^ATf ICP> 
XY(9)  *  DXK 
XY(«)  *  v?!N 


THIS  SEcTjnN  SETS  CERTAIN  VARIABLE?  FOP  ’HE  CYlINCE*  CASE. 

ASA  is  A  CONSTAT  USEC  IN  THE  SCATTERER^'S  AND  OELT*  IS  SET  To 
ZERO  FOP  *  VERTICAL  CYLINDER. 

IF(IPa  .Nf.  2)  50  TO  25 
DELTA.?. 

AKAsH  i*A</2. 

25  CONTINUE 


the  in»ut  anslfs  ah  converted  TO  PA'MANS  AND 
their  sin?*;  and  cosines  ape  calculated. 

Al PHAaAL»HA/Pir 
OELTA*r.ELTA/PAD 
CfiS0*CCStnELTA> 

SIN0«SIN(DELTa) 

COSA«rCS(ALPMA) 

SI.NA*elM(  ALPHA  ) 


BECAUSE  CP  CERTAIN  APPROXIMATIONS  -AOE  IN  The  ANALYSIS 
THERE  IS  *  LIMIT  on  the  SIZE  of  Thp  SCATTERED  That  may 
BE  SIMULATED.  TO  AVOID  THIS  PROBLEM  AS  MUCH  *S 
POSSIBLE.  FOR  THE  RECTANGULAR  SURFACE  * 

THE  PROGRAM  HILL  BREAK  UR  TOO  LARGp  A  WALL  INTO 
SMALLER  PIECES,  to  AVOID  PRORLEMS  with  other  types 
OF  SCaTTERFRS  the  VARIABLES  involveo  are  set  to  default 
VALUES  and  the  BREAKING  up  SECTION  IS  SKIPPED. 

IH*1 

iv«l 

ox*f. 

0Y*§. 

OZ»i. 

DXl»*. 

OYl«R. 

IFCIO*  ,NE.  1)  GO  TC  25 


TEMP  IS  The  MAXIMUM  DISTANCE  FROM  THE  REFERENCE  POINT  ON  ThE 
HALL  THAT  -ILL  GIVE  A  REASONABLE  ERROR  IN  THE  APPROXIMATION. 

TEMR«SORTtLAMBDA*SORT{(XXA(l)»XH(l)  )**2*(YA-XW<2n»«2)  >/5. 


I  i  L||J  JWLiai  CiiMliMTM 


noon  ciooooooo  nooonoo  n  o  oooo  oooo 


WF I TE ( ? » 3  :’13 .  T0.XW(i;,xw(2)»Xrt(3>  # ALPHA, DELTA, WW, MW.  I  m  ,  I  v 
AND  HW  ARE  SET  1  >,EW  VALUES,  T^ESE  A«E  TH£  SIZES  OF  The 
» I"CES.  "*v  AMH  OV  ARE  the  CHANGE  in  X-  ANH  V-C00«0lNATES  BETWEEN 

pieces  i*1  the  uqpihontal  tnws.  d?  is  the  change  in  elevation 
BETWEEN  *»ihs  vertically,  oyz  and  dyz  are  the  CHANGE  in  X  ano  y 

a^TWEE’'  Rf  "S.  TWJS  CHANGE  OCCURS  ->NLY  IN  TILTED  WALLS  (SIND 
NOT  EC'JAl.  TO  ZPRC). 


A  I  *  I u 
WW*Ww/ A  I 

TEHp«ww*(4l-l,>/2, 

0X*a3S  (C0SA*WW) 

Xw(l>aXw(i)-ABS(COSA*TEHP) 

0Y*SIGMSINA*WW,XW{2)) 

XW(2)aXW(?)*SlGN( (-SINA#TEMP),VW(2) ) 

HW*MH/ELCAT<IV> 

o?«coso*hh 

OXFaST-vDtHWwSlNA 
DY2«SI.ND*iW*C0SA 
GO  TO  27 


xw  is  the  coordinate  vector  used  for  the  location  of  the 

REFERENCE  POINT  OF  EACH  PIECE  OF  the  wall.  XWf*  IS  USED 
AS  ORIGIN  OF  THE  WALL.  AS  EACH  PIECE  IS  USEC  FOR  THE 
SCATTERING  XW  IS  INCREMENTED,  XW*  tS  USED  TP  RESET  Xw 
FOR  LOOPING  ON  ROWS, 

26  WRI TE ( 6 , 1013 )  10, XW(1),XW(2),XW(3), ALPHA, DELTA, WW,HW 

27  Xwm)*XW(D-DX-0X2 
XW-M2)*XH(7)-dY-0Y2 
XW3(3)*XW(3)-0Z 


this  loop  is  for  the  rows 

00  42  18*1 «  IV 
XWG a  )«xwp! (1)40X2 
XW«(2)*X*'ft(2)+0Y? 
XH£(3)*XW0(3)+DZ 
XW(i)sXW3(l) 

XW ( 2 )  aXW3 ( 2 ) 
XW<3)*XWi(3> 

C 
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-  EFn  S0UBCF  STATEMENT  .  JF'KS)  - 


C  PIECES?0  S  WI™IN  EACH  R0W  A|y°  Is  F0B  w0BI20NTallT  SFPARaTED 
•  C 

00  41  IM1.IH 
C 

C  pltcl  SlNcTSuLimX!818"  ,r  ™E  REre#EI<C'r  °M'!T  OX  THF 

c 

XW<l>sXKO  )*0X 
XU<2)sXW(?)+0Y 
C 

l  1?  Sg^Kie^r  r2ftR  SOT "  ** 

C  rwp.riEL°^  AT  ™C  REFEPENCE  b0INT  FOP  ALL  AnT£MNAE  ARE  IN 
C 

CALL  fLC(XW<i),XW<2>»XW<3)> 

c 

C  THIS  LOOP  IS  0,v  The  ANTENNAE.  FOP  PACN  PJFFF  THF  PRfiPBiM 

c  calculates  THF  scattered  field  from  all  wlLlt  0CR 

c  IFL  IS  THF  NJMPER  OF  THE  ANTENNA  BEING  SIMULATED! 

c 

c 

DO  40  IEL*1.NFL 
C 
C 

C  ANTiSna*  kR£  ™£  X”’Y"  AN°  Z”  C0M,?InatES  of  THE 
C  •  •  * 

XA  •  XXA( TEL) 

HA«1A(IEL> 

C 

\  ^HITTOEcTiOn  INITIALIZES  THE  RECEIVER  PflfNT 
Q  t0C— I0N  VAR‘Apt**$*  IP  IS  TU£  NUMBER  OF  The  receiver  point, 

IPR0 

!Ft ICr. £!),*)  go  TO  29 
COEG  *  XMIN  -  DXR 
GO  TO  30 

29  XR  a  XMIN  -  OXR 

30  CONTINUE 

IFIMODE.GT.S)  22  *  2A ( IEL ) 

C 

l  REFERENCE  ?S!J??NT‘L  D,ST*''C6  Tl(E  *"««.  TO  THE 

C 

OH  *  S0RT((XH(1)-XA)*#2  ♦  (XW(?)-YA)*#2) 

c 

\  WWia-M."  SMS‘„ 
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-MY  -  EFN  SOURCE  STATEMENT  -  IFN(S) 

c  r*r  AfTF-.-  a.  the  reference  S'-ste*  u«eo  is  alig^E'*  mTh 

C  t,je  IT c  °  •  *‘i  PECTAWGLE  *NC  TW£  THIRD  AVIS  IS 

c  The  .T'Ti»a«-  *..'OMAL.  If*  The  CASE  THE  ^YLI‘J0E»  Th£ 

C  MO^HAL  Iq  aSS-J**EH  TO  LIE  IV  A  HOPpONTAL  PLANE  AND 
C  T*  P0I'!T  AT  T-i"  A'TE'vVA. 

C 

IF <1 ~  .V.  9)  GO  TO  3’ 

AMj  )  s  ( XA-VW  ( 1 ) )  /DVc 
A\(2)s(YA-yW(2) )/CW 
ANlSIsT. 

GC  TO  33 
32  COUTIM.'E 
AN<l>sSIN* 

AN(2>s-C0->a 

A*  C  3  > . 

*3  C^vTIa:UE 


the  Hf-RIPOUTA:.  ANGLE  BEWTEEN  T*E  NORMAL  TO  THE  SURFACE  AND 
the  line  :r  sight  to  the  antenna  is  gamma,  sivg  and  cosg 

ARE  THE  ci*iE  AMO  COSINE  OF  GAMMA. 

CTiSG*.' -A v (3  )*(XH(1)-Ya)-AN{2)«(XW(?)-YA) 

SI^g  *  <-AN(2)«(YH(l)-XA)  ♦  AN(1)*(XI>,U:-YA))/0w 


IF  THF  CO^G  IS  NEGATIVE  THEN  TW£  LINE  OF  SIGHT  IS 

thru  the  rack  of  the  scatterep  and  the  illumination  of 
the  front  SURFACE  IS  ASSUMEO  tG  BE  OF  2ER0  INTENSITY 
and  THE  FILLO  from  THIS  SCATTERING  IS  IGNORES. 

IF  { Cr-Ss )  34.34,35 

34  WRITE  <6. 1017)  IA.I8.IEL 
GO  TO  40 

35  CONTINUE 
C 

c 

C  THIS  IS  The  LOOP  BACK  POINT  FOR  the  RECEIVFR  POINTS, 
c  FOR  EACH  PIECE  OF  SCATTEREP  AND  FO»  EACH  ANTENNA 

C  THE  PROGRAM  calculates  all  the  fields  at  all  the 
C  RECEIVER  POINTS  BEFORE  GOING  ON  TO  THE  NEXT  PIECE 
C  OR  ANTENNA.  XR.YR,  ANO  ZR  APE  THE  COORDINATES 
C  Of  THE  RECEIVFR  LOCATION.  VX.VY  ANO  VI  APE  THE 
C  VELOCITIES  IN  THOSE  DIRECTIONS,  THE  LOCATION 
C  IS  DETERMINED  BY  SLIGHTLY  DIFFERENT  METHODS  DEPENDING 
C  ON  THr  flight  TYPE,  the  VALUE  OF  ICF  IS  THE  CONTROL, 
g  IP  IS  THE  ^ECFIVER  point  number  ANO  is  useo  to 
C  DETERMINE  WHERE  the  FIELDS  FROM  The  SCATTERING 
C  ARE  To  BE  RUSHED. 

c 

36  CONTINUE 

IFIICF  , LE.  »>  GO  TO  37 
CDEG»CDEB*OXR 

IF (  (rOEG-XMAX)«OXR  .GE.  0.)  GO  TO  4« 
XR«R*CO$(Cr>EG/RAO) 

YR«R*SIN<CDEG/RAD) 
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4 1 N  -  £Ff.  SOUPCF  STATEMENT  -  !F\(S) 

vV  .  .  v"l»YR/P 
vy  •  V£l*XR/P 
V2  * 

GO  TO  19 
37  C0.',‘TI  ::.t 
XP*XR*CX9 

IF<  (vS-)f'MX)#DX5  .GE.  7,)  GC  tq  4P 
YP*X»*Ta‘/P= 

?b  s  »UP 

TFfXP  ,LT,  XTH)  GO  TO  Zm 
2P  *  5»9  *  <X9-XTH)*TANRR 
V?  1  VfL«SPS! 

50  TO  1,9 
39  V?  a  ?■  •  0 
39  C0«TI/'"B 

IF(IP  .GT.  490)  GC  TO  4e 
19*19-1 


-3. 

m 

I 
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PW  is  THE  ?  I  STANCE  rp0«  TW”  °FCEI VPR  PrjVT  yr.  TMF  1 

SCATTER  NEWPENCE  =01  NT.  | 

RW*SQ9T{(XR-XW{1))**2+{YR-XW{2))**?*(Z0-XH{3))*«2)  \ 

•1 

1 

RR  IS  The  HORIZONTAL  distance  from  THE  RECEIVER  to  The  I 

REFERENCE  POINT.  \ 

RR*SC07f (XB-XW(1) )e«2*CYR-XWC2 ))••?)  \ 


BETA  is  The  HORIZONTAL  ANGLE  <5rTW£FN  THF* S'lPFArE  NORMAL  and 
the  LINE  CF  SICHT  to  the  RECEIVER  POINT.  SP'R  AMO  CCS? 
are  the  sine  and  cosine  of  beta. 

C0$B*(AN<l)MXR-XW<l)>*ANm*<YR-XW(2>>>/»P 
SlNB*f-AN{2)*(XR-XW<l)  )*AM(1>#(YR-XW^) )  >/DR 


or  is  the  distance  from  the  antenna  to  t«e 

REFERENCE  POINT  ON  THE  SCATTe-PER. 

OR  a  SDRT( (XR-XA)**2  *  (YR-YA)#*2  ♦  <ZR-Z2)**?  ) 


THIS  SECT’CN  evaluates  TwF  scattering  FROM  the  SURFACE. 

THE  COMPLEX  variable  ,FAC,  RERRESEr TS  tHE  GaI*  FACTOR  3 

FROM  THE  REFERENCE  POINT  Gf)  THE  SURFACE  TO  THE 

RECEIVER  POINT.  J 

I 


PHIO  AND  PHI  Jr-  ARE  the  RELATIVE  FPFOUF”CY  SWIFTS  y.'t  TQ  DOPPLER  I  -1 

EFFECT  FROM  TuE  AIPC°AFT  VELOCITY.  %  | 

°HlO  a  A<*(VX*(XR-XA)  ♦  VY*(YR-YA)  ♦  V2* (2R-22 ) ) /OR  | 

i  ') 


’••.-*  I \  -  EpN  SOU»Cr  STATEMENT  -  JFM(S)  - 

pHljO  =  4»v*(VX«(vs<-XW(i)  )*VY*(y(?-XK(2))'t'V2#(Hi?-Xt-i(S)))/RW 


THESE  CONSTANTS  AP£  THE  GAIN  FACTORS  FOR  THF  VA&ITjS  CROSSTALK 
CASES. 

UT»(°Hl  JD-PNID)*rA/2. 

SMCUCfl)  S  SlPCi JT+a9?-T)#»? 

SMCUCf2)  s  S INC ( U7+N15 4l )**2 
SMCUT  e  SrNCUJT) 

SNCU')  «  SrNC(UT*W60T) 


this  section  calculates  the  sain  for  the  actual 

SCATTERING. 

A*AK*(SiNn#(COSG*CO$B>*COSD«<<XWC3>-HA>/OW*<XW(3>-2k>/RR)> 

3sA*?.«A<*HA*C0SD/DW 

FAC«C>:XP(CMPLX(0.,RH«AK))»UCEXP(CMPLX(i?.,A*HU))-(l.»0.))/A- 

,CEXP(CHPLX(0.,2,*AK*HA«XW<3)/O'.O>*(CEXP(CMPLX(?.,B*HW)>-(1.,R.)) 

./R> 

FAC«FaC/RW 

B*9-A 

A*AK*(SIN0*<C0SG*C0S8)*C0S0*I <XW<3)»HA)/DW»(XW(3>*?R)/RR) ) 

B*H*A 

RHP«S;PT(0P«RP*(-2R-XW(3))**2) 

FaC«FaC-(CEXP(CMpLX(0. ,RWP«AKn*{(CeXP(CHPLXCP.«A*HH)>-(l.,0. >)/A- 
.CFXP(CHPLX(0.,2.*AK#HA#XW{3)/O*'))»(CPXP(CMPLX(?,,9»HW))-(1,,0.)) 
,/B) J/RNp 

FAC»-FAC«AK»HW*C0SC/Pl/2. 

all  statements  for  calculating  the  scattering  from  rectangles  and 
cylinders  ARE  the  same  WITH  The  exception  of  THE  FOLLOWING  step. 

IDA  IS  ONE  for  the  RECTANGLE  and  TWO  FOR  THE  CYLINDER. 

IFUDa  .EQ.  i)  FAC*rAC#C0SB*SlNC(AK*WH#(IING-SlNB)/2.) 

!F<!DA  .EG.  2}  FaC*FAC»BESF<AKA,C0S8»SINR>/2. 


IF  10  IS  NEGATIVE  THE  GAIN  IS  TAKEN  IN  THE  OPPOSITE 
SENSE. 

IE<  ID  ,LT,  0)  fac*-Eac 


THE  GAIN  IS  MULTIPLIED  BY  THE  SIGNALS  AT  THE  REFERENCE 

POINT  TO  GIVE  THE  SIGNALS  AT  THE  »FCEIVE*.  THESE  SIGNALS  ARE  COMPLEX 

MAGNITUDES.  EP  IS  THE  SIDEBAND  PORTION  OF  THE  CAPRIER 

PLUS  SI0EBAN9  FOR  THE  COURSE  ANTENNA  AND  EE  TWg  SIDEBAND 

ONLY.  EM  IS  THE  SIOEBANO  PORTION  OF  THE  CARRIER  PLUS  SIDEBAND 

FOR  THE  CLEARANCE  AND  EC  THE  SIOEBANO  ONLY. 

EP  ■  FAC#EWR<IEL»i> 

EE  ■  FAC*FWR(IEL»2> 

EM  «  FAC»EWRUEL.3) 


-  EFN  SOURCE  STATEMENT  -  IFN(S)  - 


EC  a  FAC*FWR(IFI  -4) 


these  are  the  complex  phasors  for  the  signals  at  the  receiver 

POINT  FOR  THE  DIFFERENT  ANTENNAE  AND  FREQUENCIES. 

THEY  HAVE  THE  FOLLOWING  SIGNIFICANCE  I 
symbol  USAGE 

ZJp  CARRIER  FROM  THE  COURSE  ANTENNA 

2JPC(1)  90  H2  SIDEBAND  FOR  COURSE 
2JPC(2)  150  HZ  SIDEBAND  FOR  COURSE 
|JM  CARRIER  FROM  CLEARANCE 

|JhC(1)  9P  H2  FROM  CLEARANCE 
2JMC(2)  1?P  HZ  FROM  CLEARANCE 

ZjP  a  EP/CMPLX(CP.e.«) 

ZjPC<l)  a  EP  -  EE 

2jPC<2>  a  EP  ♦  £E 

2jM  a  EM/CMPLX(CM,0.0) 

ZJMC<1)  a  EM  -  EC 

2jMC(2)  ■  CM  ♦  EC 


SUBROUTINE  VABCAL  ADOS  THE  FIELDS  Tft  THE  FIELOS 

accumulated  pcr  the  .IR.Th  receiver  point, 
call  varcal  hr> 


THE  PROCRaM  LOOPS  BACK  TO  thE  NEXT  RECEIVEP  POINT. 

GO  TO  36 

40  CONTINUE 

41  CONTINUE 

42  CONTINUE 


THIS  IS  The  TRANSFER  BACK  TO  PICK  UP  THE 
NEXT  SCATTERCR  OR  CONTROL  CARD. 

GO  TO  21 


at  this  point  the  program  has  accumulated  the  scattered  fields 

AND  HAS  RFAO  in  A  CONTROL  CA»0  TERMINATING  THE  RUN, 

THE  PROGRAM  WILL  ADD  IN  THE  01»ECT  UNSCATTERED  FIELD.  BOTH 
DIRECTLY  FROM  THF  ANTENNA  AND  REFLECTED  FROM  THE  GROUND, 

THEN  .THE  APPROPRIATE  RECORDS  KILL  «E  OUTPUT. 

43  CONTINUE 

IR*B 

SNCUT  a  1.0 
SNCUD  a  0.0 
SNCUC(l)  a  «. 

SNCUCC2)  a  «. 


. *'*■*'“ . . . . . * . * 


1  \ 


£FN  <?CuRrr  STATEMENT 


if\ks) 


04/J 


C  FROM  TMJS  r* 7  a  TFMFnT  THROUGH  JiJST  arr^SF  ^TaTFMFNT  pi  IS 
C  THE  unop  o:--  RECEIVER  POINT.  TH£  L*0»It.G  IS  Of>E  T|-t  SAME 
C  AS  THF  SfTTi:-,  FOLLOWING  STATEMENT  3S. 

c 

44  I F ( I Cff  .ST.  0)  GO  TO  46 

xr  ■  yf  ♦  •■'xp 

I F  (  ( Y'-y  *AX)«f)yP  .r,E.  £* .  )  GO  T7  fj 

YS  s  V'«T<,:.PR 

IF(XN  .LT.  XTh )  GO  Tr  ;5 

2R  a  (  xR-vtH)*TANS!>  *  -«ua 

V2  b  vel*spsi 

GO  TO  47 

45  V*  b  3. 7) 

*R  b  ?L'P 
GO  TO  **7 

46  COEG  *  COES  ♦  OXR 

IF  <  (CPF G~XMAX)*DXR  .GE.  7i.  )  GO  TO  51 
TE^P  s  OJLu/PAO 
XP  b  p*C0S(TEMP) 

YR  b  P»Sr:(TEHP) 

47  IR*IR*1 

CALL  '•LEA°(2E,4) 

C 

C 

C  THIS  CALL  TO  flc  CAUSES  THE  CALCULATION  OF  THE  FIELD  LEVELS 
C  AT  THr  RECEIVER  POINT, 

C 

CALL  FLCOR.YP.2R) 

C 

C 

C  THIS  IS  ThE  LOOP  FOR  THE  OIFpEPENT  ANTFNNAE.  IEL  IS  ThE 
C  ANTENNA  NUMBER.  NEL  IS  TOTAL  NUMBER  OF  ANTE^naE  BEING 
C  USEO. 

C 

00  49  IElbi.NEL 
C 

c  this  section  calculates  The  fields  fop  THE  VARIOUS  signals 

C  AT  THE  RECEIVER  POINT. 

c 

HA  b  2 A { IEL) 

XA  b  VXAIIFL) 

P0»SQ»T(RaUEL)*<*2-(2P-HA)**2) 

CE*CMPLX(PO/PA(IEL),e.) 

RD*2.*AK»MA«2R/R0 
CE«CE*CMPlX(1,-COS(RO),-SIN(R0) ) 

DO  5?  J  =  1.4 
EWRdFL.  j>bEWR(IEL.J)*CE 
50  ZE(J>s2E(J)*EWRUEL.J) 

HjP  •  EWR' IEL.l)/CMPLX(Ce.2.0) 

2jPC<<>  *  EWRUEL.l)  -  EUR ( IEL. 2) 

2jPC<?)  «  EWR{  IEL.il  ♦  EWRUEL.2) 

2jM  b  EwPCIEL,3)/CMPLX(CH.0.i) 

2JMCID  *  EWR( IEL.3J  -  EWP( IEL.4 ) 

2jMC<?>  ■  EWR(  IEL. 3)  ♦  EWPUEL.4) 

e 

C 


42 
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annoonnnnnnn  ooooo  ooooo  oooooo  onnononn 


9)4/1 

MAIN  -  EFN  SOURCE  STATEMENT  -  JEMS)  - 

C  *HIS  CALL  TO  vARCaL  aoos  the  HELDS  TO  T*e  ONES  ACCUMULATED 
C  FROM  THE  SCATTERERS. 

c 

CALL  VARCAL  (IR) 

49  CONTINUE 


detec  takes  the  complex  fielo  phasors  and  evaluates 
the  course  DEVIATION  INDICATION  (con.  IR  is  the  receiver  point 
NUMSER  AND  is  USEO  IN  THE  SURROUTImE  TO  SELECT  WHICH  FIELDS 
ARE  TO  IE  USED.  DF(IR)  IS  THE  LOCATION  IN  T«E  ARRAY  WHERE 
THE  COI  IS  TO  PE  PLACED. 

CALL  OETEC  (IR,DF(IR)) 

IF(IR  . 6T .  499)  GO  TO  51 
CO  TO  44 
51  CONTINUE 

XY<1S)«FL0AT(IR) 

WRITER, 1*18)  ID.NC.IR.ICF 


THIS  SECTION  OUTPUTS  THE  COI  ON  UNIT  8.  THE  OUTPUT  IS  A  LABEL 
RECORO  (HFMO),  TWO  RECOROS  OF  FLIGHT  AKD  ANTSnwA  DESCRIPTION, 
AND  THE  COI  RECORDS. 

IF(I0  .EO.  0)  HEM0(13)«ILBL(5) 

HPITE  (6,1005)  MEMO 
HPITE(8, 1814)  XY.IO.NC.ICF 
WRITE f 8, 1016)  (DF(I).I«1,IP) 

IF  THE  10.  IS  NOT  0  THE  FLOW  IS  TO  STATEMENT  *7  TO  PROCESS 

the  id  value  from  the  control  card,  otherwise  the  signal 
strengths  are  output. 

IF (  ID  .NE.  i  )  GO  TO  57 


ix  is  the  number  of  signal  types  that  are  to  be  output,  two 

FOR  SIMPLE  ANTENNA  SYSTEMS,  FOUR  FOR  CAPTURE  EFFECT, 

IXM 

!F(NEL  .E5.  1)  !Xc2 


THESE  LOOPS  CALCULATE  THE  SIGNAL  STRENGTHS.  THE  VALUES  ARr 
PLACED  IN  XXRv(I,j).  WHERE  I  IS  The  RECEIVER  POINT  NUMBER  AND 
J  HAS  THE  FOLLOWING  USAGE) 

J  USAGE 

1  CARRIER  LEVEL  FOR  COURSE  ANTENNA 

2  sioebano  level  FOR  COURSE  anrenna 

3  CARRIER  LEVEL  FOR  CLEARANCE 

4  SIOEBANO  LEVEL  FOR  CLEARANCE 

XXRY  OCCUPIES  THE  SAME  LOCATION  IN  CORE  AS  ZP  AND  ZH. 

DO  92  1*1. IR 

52  XXRY( I f 1) sCABS (2P( I ) >*8.2 


o  o  o  o  o  Cl  o  r»  oo  o  oo  ooo  nooooo 


*AI‘*  -  EFn  SOURCE  STATEMENT  •  XFS(S>  - 

OC  33  1*1. I» 

53  XXRYCT,2>*CARS<2PC(I.l>-2PC<I,2))/2. 

00  54  1*1,  Ifi 

54  XYPY<!,3>*Cm<lMm>#0.2 
00  55  1*1, IR 

55  XXRv(i,4>*CA8S(1MC(I,1>-2MCU,2))/2. 


this  loop  outputs  the  appropriate  number  or  SIGNALS  ON  UNIT  8. 
the  label  record  tor  each  case  is  altered  slightly  as  explainer 
in  the  data  STATEMENT  for  ilbl.  1  * 


00  54  J  *  1,IX 
MEM0(13)*ILBL(J) 

WRITE(B,1*«5)  MEMO 
WRITE(8.1014)  XY, ID*NC, IGF 
56  WRITE(S,1016>  CXXRYtl . J) , 1*1. JR) 


this  section  controls  the  flow  of  the  proopam  after  the  ouput 
for  The  C*SE  is  finished.  THE  CONTROL  IS  BY  TmE  VALUE  OF  THE 
ID  READ  IN  ON  THE  LAST  CONTROL  CARO.  THIS  ABSOLUTE 
VALUE  OF  ID  IS  IN  IDA,  DEPENDING  ON  T«E  VALUE  OF  IDA  THE 
PRCGRaH  LOOPS  BACK  AND  READS  IN  THE  NEXT  DATA  CARD  FQR  THE 

nextcaSe  to  BE  run.  the  value  hill  cause  The  transfer  in 


the  following: 
IDa 
0 

3-lB 
11-15 
16-2B 
21-5] 


NEXT  TYPE  OF  card  RCAO 
MODE 

SCATTERER 

LABEL 

MODE 

COURSE  *IOTH 


57  CONTINUE 

IPUOA  .EO.  0)  GO  TO  l 
IFCIOA  ,LE.  10)  GO  TO  20 
IF < IDA  ,LE.  15)  GO  TO  14 
IFIIOa  .LE.  20)  GO  TO  1 
IFCIPa  .LE.  50)  GO  TO  6 

58  CONTINUE 

end  file  8 

PEW  I  NO  • 

STOP 


FT, 2/ 

3F9.2/ 

FT.2.8H  DECREES 


1000  FORMAT  (8F10.3) 

1001  FoPMAT(I2,2X,6X,7F10.3> 

1002  FoRMAT(5X#4HCLS*,F9.4) 

1003  FORMATCSH0HOOE  •  2U/10H  FRQ  a 

1  «H  XTH  s  F9.2/  8H  2A  • 

2  8H  Xa  a  3F9.2/14H  COURSE  WIDTH 

1004  FORMAT  (3X,13A6,A2) 

1005  FORMAT  (1JA6.A2) 

1006  P0RMA1  (7P10.0.2X.3I2) 

1007  FORMATtiHS  VEL«.EU,4) 

1008  F0RMATI26H  OVER  300  RECEIVER  POINTS  ) 

1009  PORMAT<6H0XM!N«,E11.4,7H  XMAX*.EllV4.7H  DXR.,|U.4,7H 

X.4«7H  PSIR*,E11.4.6H  XTH*. Ell, 4, 5H  IUP«,|11,4) 


PHIR«.E11 
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M  A  I  ** 


-  EFn 


csurcf  statement 


-  IF-'H?)  - 


ALaHA 


1819  FORMAT  (li-<r  STRUCTI'PE  data) 

1011  FORMAT C56H  10  XK  VH  2N  iLaHA 

X,0MQEl7A  ,5X,23H  KW  ‘•w  .fX.JHH 

X.5X.13H  *  SECTION?  ) 

1012  FORMAT  (12, 3e  .0.5X.2F5.F.3F10.0) 

1013  FORMAT  <I?,l<.7Ei2.*»5X,!3.4X.t3> 

1014  FORMATtiV^rif..?,/  3Fl3.«.IlO.liy,?I10J 

1015  F0RMA7-(6NcMlNn*.-11.4,7H  MAXDb,E1J ,4»7h  'iDE&B,ru.4,4M  Pa, 

XE11.4.7H  2UP*,E11.4,7H  ICFa,12> 

1016  Format (  7F15.8  ) 

1017  FORMA" (27*-  SURFACE  IS  NOT  ILLUMINATE"  . 

X5H  h*,I2,5H  Va,I2,6M  IEL*.I2) 

1010  FORMAT  (?X,3MI0*  ,I3»5X#3WNC«  .H.fX.SHlRa  , !3.5X,AHlCFa  ,12,//) 
EnO 


'iOE6«»FH,4,4M  Pa, 


4S 


ct  n  oo  o 


SL'?1 


ZFN  SOURCE  STATEMENT 


ienis) 


this  sue i ‘.*e  is  used  to 

VARIOUS  M»TR4riES. 


ZERO  OUT  THE  CONTENT*  OF 


Si.'RROuTive  CLEAR  (X,N) 

complex  x  f i > 

00  1  I  «  t,N 

i  x<n  * 

RfTURv 

E?iO 


?4/ 


SUB? 


-  EFn  cOllRCf  CT  ATEMENT  -  IF‘.*(S)  - 


c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

e 


c 

c 

c 

c 

c 

c 

c 

e 

c 

c 

c 

c 

c 

c 


THIS  SUp«OUT!ve  IS  USED  TO  T  se»JT  Data  ro«  CALCULATING  THFPPECT  I  CAL 

patterns  ros  array  type  ante^-ae. 

Sl'PROuTl*?  CRRNTS (  O*  C.  S»  ET.  *•'£  ) 

LOGICAL  EOF 
OIMEN*IO‘i  ET(i9),D(l) 

COMPLEX  CCl).S(l) 

COMMON  /SUB/  MODE.ICP.FRfJ,LAM8r.A#PI,»Ano.PwI(3).PSl{?),NEL.XTrt 
IFtEOFU))  G«  TO  3 
I  «  1 

this  is  the  input  for  the  element  location  and  current  description 
OT  IS  The  element  DISPLACEMENT  in  THE  Y-PIRECTiON.  measured 

in  wavelengths. 

CT  IS  The  CARRIER  Plus  SIDEBAND  AMPLITUDE,  in  PELATIVE  units 
PC  IS  The  ca»«ier  PLUS  sideband  phase,  in  degrees 

ST  is  The  SIDEBAND  ONLY  AMPLITUDE,  I?'  RELATIVE  UNITS 
PS  IS  THE  SIDEPANO  ONLY  PHASE,  IN  DEGREES 

1  READ  (5.10P0)  DT,  CT.  PC.  ST,  PS 

THIS  test  IS  ID  See  if  the  ENO  of  THE  ELEMENT  CAROS  MAS  BEEN 
RFACWE3.  Ir  THE  CARRIER  PM/SE  IS  GREATER  THAN  500  FLOh 
is  to  the  elemfnt  pattern  section. 

IF(  Pc  .ST.  500.)  GO  to  ? 

THIS  IS  The  #0  DEGREE  PHASE  SHIFT  FOR  THE  OUAOPATURE  OF 

the  sidepano  only  to  the  sideband  in  tme  car°ifR  plus  sideband. 

PS  ■  »s«91.l 

HRITE  (0.1000)  OT.CT.PC.ST.PS 
0(1)  *  DT*2.«PI 

C(I)»CT»CEXP(CMPLX(0..PC«RADD) > 

SCI)  =  ST<*CEXP(CMPLX(P.,PS*RAC*')) 

I  *  I  ♦  1 


THIS  STATEMENT  LOOPS  BACK  COR  THE  *>EXT  ELEMENT  IF  T«E  TOTAL 
NUMBER  OF  ELEMENTS  OOES  NOT  EXCEED  THE  AVAILABLE  SPACE. 

I F (  I  .LT.  26)  GO  TO  1 


THIS  SECT»CN  RPAPS  IN  THE  PATTERN  F0»  ThE  ELMERS,  N£  IS  T«E 
NUMBER  OF  FlEMFNTS.  ALL  ELEMENTS  aRF  ASSUME"  TO  HAVE  THE  SAME 
PATTERNS. 

2  NE  ■  I  -  1 

ET  KILL  COMT*I*t  THE  ELEMENT  PATTERN.  THE  VALUES  ARC  IN 
RELATIVE  AMPLITUDES.  ETC 1>  IS  THE  VALUE  AT  *ERO  0E6EES  AND 
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. . 


-  EPS  SOURCE  STATEMENT  .  jp^S) 


C  SUCCESS  I  V“  VALUES  ARr  *T  i**  rtrfiPFr  CDif<iui>  . 
C 

READ  (5.1?, 28)  et 
KPITE-'-StC^)  ET 
?fTU°V' 

3  W5 1  TE ( 6, 1  ?<•.! ) 

E\M?  EtlE  S 
STOP 

ippc  fcrmatc  ini.4  ) 

3!WT  *****  »if$smc  j 


I.u,  ^■LHtliL...rt„Ji  I.h  n  '.-.J!  M.  '-HM-H.I  .LKilii  l!i.i  >■!.  ■■■■.iCj:  .ly-jULl!  U.  UlJili!.!liJ..irjill.Jl.iuH.il...J..J!  lull  a,  Juto  3£MUAL4*£u/Al«.  JA  X.'Jjil  AJ.  m"wJU,!tl.,PillB,llliill!-UiL1j  jl  itl'OliluilllUJl-Hluii.li  Jililill  JIaiIiJI  ,1  JOIlILI  Jlililill  Alta  t  j  UiL  iL'kJ,  P!ii  »..JL  H,li..«.l'U.1,jLl  t'..ti,il',.llL..  liL!la..L.ii.jiim!Ii 


nnnnn 


SUB3 


EFN  SOURCE  STATEMENT 


Z*/ 


T«!|  SUBROUTINE  INPUTS  the  ANTENNA  PATTERNS  FOR  The  hEaSUREO 
P»TTERN  antenna  Cases, 

SUSROjTive  PATTRN  (  ARAO,  aEPP,  aGPP  I 
LOGIC' l  £Or 

DIMENSION  ARAD(5D ,  AEPPC52).  AS*P<5P) 

DATA  RAD  /  57.2957795  / 

IX  ■  1 

ir(E0F(5>)  GO  TO  4 

1  REA0I5.17P3)  AN5.  AFPpUX).  AC»P(IX) 

APPP< lX)«APPP(lX)«100RCe. 

ACRPClX)«AGPP(IXJ»1836>f . 

arad<»x)«anc  /rad 

IX«IX#1 

IF<  IX  .3E.  51)  GO  TO  2 
IFC  AnG  .LT.  361.)  60  TO  1 
IF C  lx  .LE.  2)  CO  TO  4 

2  WRITE  (6*1001) 

WRITE  (6*1*02) 
lY«I*-2 
DO  3  !*1» I Y 

ANGsARAOC I )*RaD*. 00001 

3  WRITE  (6.1003)  ANC»AEpP( I ) > ACPP( I ) 

GO  TO  5 

5  RETURN 

4  WRITE  (6.1004) 

end  pile  r 

STOP 

1RR0  PORMAT(0PiP.0) 

1911  PORHAT(2BH0ANTENNA  PATTERN  MEASUREMENT) 

1R02  P0RMATCJ4H  ancle  READ  SIDEBAND  CARRIER) 

1093  FORMAT  (3E12.4) 

1RB4  F0RHATC33H  MEASURED  ANTENNA  PATTERN  MJSSINC  ) 

END 
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;up< 


-  EEN  SOURCf  STATEMENT  -  I F N < S > 


C 

c 

C  THIS  SUpRiUTIN'E  SIMULATES  the  BEHAVIOR  op  the  ils  receiver 
c  SVST-H.  c"OR  THE  IR,  TH  RECEIVER  POINT  IT  CALCULATES  THE  CD  I 
C  THAT  wDULn  BE  OBSERVED  WITH  ThE  FIELD  LEVELS  IN  ZP.ZM 
C  ZPl  AND  ?!:C, 

U 

Subroutine  detec  ur.coi) 

DOUBLE  PRECISION  G000 

real  n 

COMPLEX  HP(500>,ZPC<500,2), 

2  ZMI500) ,ZMC(500.2> 

DIMENSION  VCD ( 500 » 2 ) » VPD <500* 2 ) »  VMD (300 » 2> 

DT'IENSIOM  V(2)  »GC00(26> 

COMMO'.  2®, 2PC, ZM. ZMC»  VCD. VPD. Vmq 

COMMON  /VAR/  SM,SNCUT,SNCU0»SNCUC(2>*VPC(2)»VMC(2) 

DaTa  \'G  /5/ 

DATA  G30*/  .10000O*1»-»2500000Q*00, 

1- ,  4  6  875000P0000  00-01,  -  .}.  9531 25000000*0-0 1 » -  •  106811523437500-01 » 
l-.672®1259765625D-02,-. 462627410808670-02, ".337529182434030-02. 
1-. 2571 02306932210-02. -.202349037833310-02. -.163396848074630-02. 
1-. 134701120623500-02. -.112952502189500-02. -.960764626611880-03. 
1-. 827188932350790-03.-. 719854371145180-03. -.63180593716750D-03. 
1-. 559H546169754D-03. -.498205770262850-03, -.446369856295290-03, 
1-. 4030120751646310-03. -.365323232359670-03, -.332678129480020-03, 
1- . 304221257334890-03 . - .279265607319140-03 , - . 257259477462390-03/ 
CALL  OTC(ZPUR),V?,VPC) 

CALL  OTC<0H(IR).VH,VMC) 

BK2  *  4.0*VP«VM/( VP*VM)»*2 
UW  ■  VM/VP 

IF <  UN  .EQ.  0.)  GO  TO  2 
N  *  NG 
Nl  *  N  t  1, 

CC  «  0.0 
CP  *  0.0 
CM  ■  0.0 

1  IE<  Ni  «LE,  0  )  GO  TO  3 
G  >  G000(N1) 

CC  *  CC*BK2  ♦  G 
CP  a  CP«0K2 
CM  a  CM»BK2 
Nl  a  Nl  -  1 
NaNl-l 
GO  TO  1 

2  CC»  1.0 
CP  a  1  .0 
CM  b  *,t 

3  00  4  I  a  1,2 

V02I  8  CPBCP«VPDUP,U  ♦  CM*CM«VMDUP,  n  ♦  CC#CC*VCOUR.I) 

VCI  ■  CP4VPCU)  ♦  CMBVMCU) 

4  V(I>  t  SORT (  VCIBVCI  +  VOS  I  ) 

COI  »  SM»CVC*>-V(l))/{Vt2)*V<i)) 

RETURN 
END 
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G*(l.*N#(UW-l. ) ) 

6*  <  1 .  *N*  ( 1 ,  /UW-.l , ) ) 


,4,jvk  LU-viJJJ  .^W 


kj+J--;1  iwuiggtf 


I 


r"!j?wi;i«»'tn'.- » 


sUB* 


-  EF'.|  source  S  Ta TPMENT  -  X P 'J c «5 ) 


c 

sjvjms" 

SUBR0.;Tp;r  0Tc  (  £N,VN,  VNC  ) 

cnMPtrv 

DIMENSION  2*1(500. 1),VNC<1) 

VN  a  CABS(?N(1,1) > 

Pm  a  ■*,{» 

cSspVl:  inI(*pS;>  BH  *  AT*^'2<  AImaC(2N'(1,d  )  ,BEAL(iN(l#  i  > ) ) 

<?INPa<5iN(oU) 

W  ]  I  e  «,2 

1  RETURN  S  C0SP#R£AL(ZK,‘1'I  +  U>  ♦  SlMpaMMAC(2N(i,i*i)) 
end 


51 
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SU96  -  £FN  SOURCE  STATEMENT  -  IFN(3>  - 


THIS  pU8R0UTJNE  AODS  THE  FIELDS  IN  ZJP,  IJM,  ZJPC.  AND  2JMC 
TP  THE  SUMMATIONS  IN  2*0.  HMC,  VCD.  VPO  AND  VMD*  THE  ARRAYS 
CONTAIN  TmE  COMPLEX  SUMS  FOR  EACH  RECEIVER  POINT,  ThE  SYMBOLS 
HAVE  THE  FOLLOWING  USAGES 
SYMBOL  USAGE 

ZP  CARRIER  FROM  COURSE  ANTENNA 

2M  CARRIER  FROM  CLEARANCE 

ZPCUR.l)  90  HZ  SIDEBAND  FROM  COURSE 

IPrJ  IR.2)  150  HZ  SIDEBAND  PROM  COURSE 

ZMC(IR.I)  90  HZ  SIDEBAND  FROM  CLEARANCE 

ZMCUR.2)  150  HZ  SIDEBAND  FROM  CLEARANCE 

VCDUR.l)  • 

VC0UR.2)  * 

VPOdR.l)  *  THESE  ARE  INTERNAL  VARIABLES  USED  FOR 

VPDMR.2)  •  DOPPLER  EFFECTS,  THEY  HAVE  NO  DIRECT 

VMO(lR.l)  •  PHYSICAL  MEANING. 

VMq (IR#  2)  * 

SNCUT  IS  THE  GAIN  FACTOR  FROM  THE  DIFFERENCE  OF  THE  SCATTERED 
SIGNAL  FROM  THE  DIRECT  SIGNAL  FREOUENCY.  THIS  FREQUENCY 
SHIFT  IS  CAUSEO  BY  THE  DIFFERENT  VELOCITIES  OF  THE  AIRCRAFT 
RELATIVE  TO  THE  ILS  ANTENNA  ANO  the  SCATTERCRS.  SNCUCll)  IS 
THE  GAIN  OF  THE  CROSS  TALK  FROM  THE  CARRIER  THROUGH  THE  90  HZ 
FILTER.  SINCUC<2)  IS  THE  CROSS  TALK  AT  15R  HI. 

SNCUO  IS  THE  CROSS  TALK  FACTOR  BETWEEN  THE  90  Hi  AND  150  HZ 
SIGNAL*  FROM  THE  DOPPLER  SHIFT. 

SUBROUTINE  VARCAL  (IR) 

COMPLEX  Z 

COMPLEX  ZP(500).ZPC(5B0i2). 

2  ZM(500) ,ZMC (500.2) 

DIMENSION  VCD ( 500 . 2 ) . VPO ( 500 * 2 ) . VMO ( 090 . 2 ) 

COMMON  ZP.ZPC.ZM.ZMC.VCD.VPO.VMO 
COMMON  /VAR/  SM.SNCUT. SNCUO. SNCUCC2) 

COMPLEX  ZjM.ZjP, ZJPC(2) »ZJMC(2) 

COMMON  /«B/  ZJM.ZJP, IjPC.ZJMC 
CA02<Z>  •  REAL (I*CONJG(Z) ) 

ZP(IR)  ■  ZP(IR)  ♦  ZJP 
ZMUR)  •  ZM(IR)  ♦  ZJM 
00  1  1*1.  2 

ipcdR.i)  •  zpedR.n  *  zjPcd)*SNcur 
1HCUR.I)  •  iMCdR.I)  ♦  ZjMCd)*3NCUT 

VCO(tR.I)  •  VCOdR.D  ♦  (CAB2(ZJP  )  *  CAB2(ljM  )5«SNCUC(I) 
J*S-I 

SNCUO*  •  SNCUD«SNCUO 

VPOUn.J)  •  VPDdR.J)  ♦  CAB2(lJPC(  I ) )  «8NCU02 
1  VMO(IH.J)  ■  VMO(IR.J)  ♦  CAlS(ZjMCd))  «SNCU02 
RETURN 
END 


S2 


»»',i  i  *v  im  v* 
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-  EFN  SOURCE  STATEMENT  •  jFn(S)  - 


this  su8»*5Utine  calculates  the  electric  fields  pob  the 
sidebands  AT  LOCATION  (XI, Y. 2).  array  e  is  the  same  as 
array  ew»  in  the  main  PROGRAM. 

subroutine  flc(xi.y.z> 

COMPLEX  £.P,FPP.GPP.C(25.2>.S(25.2> 

COMMOrJ/CD/  ARAO<53),AFPP(50>.AGPPt50>*RRAD(50>.0FPP<50>,8GPP<5F> 
COMMO-.J  /5u«/  LC(2>.FRQ,NAN0A#P!.RAnDiPHI.P{2),P$I,TT<2>*NEL»XTM, 

1  XXA(3)*YA.HA(3).RA(3) 

COMMON  /ANT/  L0C.FPP<2),GPP<2>,E(4,4>»CHA(2>.AS(2>,D<25»2>.C.S, 

2  ET(20,2) »ND(2) 

AK*?.*PI/HAMDA 

JA  *  i 

this  is  the  loop  on  antenna  number. 

DO  1  J*l. NEL 
CALL  "LEAP  CFPP.4) 


LOC  IS  THr  ••'PE  FOR  ANTENNA  .J. 

LOC  ■  LC(J) 

X  IS  the  distance  FROM  the  antenna  TO  THE  POINT. 

X  •  *t  -  XXA { J  > 

R«SQRT<X»«2+Y#«2*<2-HAtJ>>»*2> 

RA(J>  *  R 
PHl«ATAN2(Y,X>. 

PSI  •  ATAN2(I-HA(J)#X) 

JA*l*l#f*(J-l) 

IF(  LOC  ,LT.  4)  CALL  CSP 

IF(LOC  .CO.  5)  CALL  ANTP(FPP( J) .GPPC J) . ARAO< JA ) , AFPpt JA) , AGPP( JA > > 
IF (LOC  .EQ.  7)  CALL  LNAR(FPP( J) .GPP( J) .PHI .0(1. J) .C(l. J) . 
.S(1,J).ET(1.J).N0(J)) 

CONS  «  AK»R 

F  IS  THE  COMPLEX  GAIN  FACTOR  FOR  THE  TRANSMISSION  LOSS  FROM  THE 
ANTENNA  TO  THE  POINT. 

F  •  CEXP(CMPLX(f ..C0N3U/R 

00  1  jCal.2 

JBrR«jC-l 

GPP  IS  the  signal  LEVEL  For  the  SIDEBAND  PORTION  of  THE  carrier 
PLUS  SIDEBAND. 

GPP(JC)*  GPP(JC)«AS(JC> 


law  iiy*  hi*  w  W  irw  kWiSjjl? 


|P  - 


► 


su>*  -  EFN  SOURCE  STATEMENT  -  IFN(S)  - 

C 

C  FPP  IS  THf  COMPLEX  PHaSOR  FOP  THE  SIDEBAND  ONLY, 

c 

FPP(JC>«FPP( JC)*CHA( JC)*AS<JC) 

E ( J » J9  > «6PP  <  JC ) »F 
1  E(J,JB+i)s  FPP(JCJ«F 

return 

END 


n 


*IDPTr  SURP 


TH!<s  At IRROIJT I NE  GIVES  FPP  AND  GPP  AT  ANGLE  PHI  BY  SUMMING  THE  SIGNAL: 
C  FROM  THE  ND  ELEMENTS  IN  THE  ARRAY.  THE  PATTERN  FOR  THE 
C  FlEMENTS  IS  IN  FT.  THE  RELATIVE  CARRIER  PLUS  SIDEBANDS  AND 
C  SIDFBAND  ONLY  SIGNALS  FED  TO  THE  ELEMENTS  ARE  IN  C  AND  S. 

C 

SUBROUTINE  LNAR  ( FPP ,GPP .PHI »D »C ,S »ET »ND ) 

COMPLEX  FPP.GPP.C.S 
DIMENSION  D( 1) »C( 1) ,S(1) »ET(IJ 
S I PH=S IN { PHI ) 

T EMP=ABS< PHI )/. 1745329 
I=TEMP+1. 

A  =  I-1 

P=TFMP-a 

EPP*R*(FT(  1  +  1  )-ET<  in+ETd) 

FPP= (O.O.O.rt) 

GPP=(0. 0,0.0) 

DO  1  J*1 ,ND 

GPP  =  GPP  +  C( J)*CEXP(CMPLX(0.,-D( J)*SIPH) ) 

1  FPP  *  FPP  +  S( J)*CEXP(CMPLX(0.»-D( J)*SIPM) ) 

GPP  *  EPP*GPP 
FPP  *  EPP*FPP 
RFTUPN 
FND 
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oooooooo 


'Uao  -  £F‘.  SOURCE  STATEMENT  -  ir-g(S) 


*4/ 


THIS  Af'TFNNA  SUBROUTINE  GIVES  PPP  AND  GPP  FOR  ANGLE  ®HI  SY 

INTERPOLATION  in  tables  ant  and  acp.  angle  phi  is  in 

RADIANS.  THE  SUBROUTINE  WILL  INTERPOLATE  BETWEEN  VALUES 

pbacketthg  tphi.  if  phi  is  outside  the  Range  of  the  ta»lf 
then  extrapolation  from  the  last  two  VALUES  will  PE  used. 

SUBROUTINE  ANTP  (FPB.GPP.ANG.ANT.ACP) 

DIMENSION  ANGC50) .  ANT($&>,  ACP(3B> 

COMMON  /SUB/  LC(2)  .FRO.WAMDA.pi .RADOiPHI ,P(2),0Si.T(2),NAR,XTh, 

1  XXA(3),VA,HA(3).RA(3) 

D n  l  1*2. 5C! 

Ka  I 

IF(ANG(I)  ,CE.  &.S)  GO  TO  5 
IFUN.^m-PHI)  1.3,2 

1  CONTINUE 

2  FPP«A*jT(*<-l)*(ANT(K)-ANT(K-l)  )»(PHI  •AMSfK-l)  )/(ANG(K)-ANG(K-i>) 
GPP«ACP(K-1)*(ACP(K)-ACP(K-1))#(PHI  •ANG(K-1))/<ANG(K)-ANG(K-D) 
GO  TO  4- 

3  FPP»ANT(<> 

GPP«ACP(K) 

4  RETURN 

5  K*K-1 
GO  TO  2 
END 
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onno  anodao 


'UB10 


-  efm  sourer  statement  -  i r ^ c s >  - 


*4/; 


THIS  ANTElA  SUBROUTINE  WILL  EVALUATE  FPP  ANQ  GPP  FO*  THE 

standard  antennae.  the  value  of  lcc  will  oeterhine  THE  type 
OF  ANTENNA  USEO.  the  SIGNALS  WILL  re  calculated  at  angle  phi, 

Sur‘P0*'TlNC  CSP 
Real  u ah3a 

CCHHOv  /SUB/  LC(2>#FRQ..WAmDA,PI,RAD0iPHI.P<2)»P$I»TC2>.NAR,xT* 
COMMON  /ANT/  lOC,FPP»XF,FPM,YF,GPP»X6»GPM,y6»£(4.4) 

dimension  c(i0).sa0>,o<ie>iET(2«p> 

SlPH*SINtP«I) 

GO  To  (1 . 4#6) *L0C 


THIS  is  The  v-ring  antenna 


Cf*2.22i 

C(1)«1.02B 

C(2)*0.546 

C(3)«(*.3S5 

C(4)«pi.2?5 

C(5)*p.214 

C(6)«?.175 

C(7)*3.i48 

0(1>*186.9 

D(2)*497.4 

0(3>*786.8 

0(4) *112?, 

0(5)«1442, 

0(6)*?.763. 

D(7)*2C|3. 

DO  2  J*l»7 
D( J)*0( J>»RADD ’ 
ET(1)*1.8« 

ET(2>*P,99 

ET(3)*».97 

ET(4)*f.f2 

ETCS)*?, 14 

ET(A>s0.7M 

ET(7)*fi,62 

ET(8)*0.4a 

ET(9)*0.S3 

ET(l2)*|,22 

ET(ll)*i,13 

ET(i2)*|,l3 

ET(l3)«f ,18 

ET(14>*8.23 

ET(H)«S,3P 

ET(l8)*f,36 

ET(l7)««,38 

ET(1B)*8.39 

ET<1*>M.40 

TENP»ARS(PHn/. 1745329 
!«TEHP*1. 
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SUB10 


EF^<  SOURCE  STATEMENT 


IFN(S) 


C 

C 

C 

C 


A  =  I-1 
R«TEMP-4 

EPP«R»(ET{l4i)-ET(i))4rTU) 

rp Pa?.? 

GPP«C?»EPP 
DO  3  J*l»7 
CSPM*nOS<D<J)#$IPM) 
SNPM*S!N( V(J)#SIPH) 


-  Ho  ■  ill  *  •.-E-WJI.CSWI 
?  I!  .!  ♦  2-#EpP#C(J)*SNRH 

50  TO  A 


THIS  IS  T^E  8-LOOP  antenna 


c 

c 

e 


4  C(l)*l«19 
C(2)*i.«| 
C(3)ap.M 

C(4)«0.3J 

0(l)«55.l 
0(2)«i9f .• 
0(3)«5*I.P 
D(4)«9Cf.(| 
FM«I. 


CSPHp2.«COSCRaOO»0(1)»SIPh> 

GM«C(1)«CSPH 

00  5  J«2#4 

I N  (  RaOO«C  < J  )  #S  I  PH  > 
FpP«rpP*C(J)*SNRH 


GO  TO  • 


c  this  i*  the  waveguide  antenna 

c 


6  C(l)«3,2t> 
C(2)P2.9fp 
C(3)«2.Stf 

C(4)»2.|2| 

C(5)«l.4l* 

C(i)*t.uj 

G(7)«i,f 45 

C{I)>MM 

C(9}«.R.l6 

Sm*|.179 

S(2)a«.913 

S(3)M.774 

S(4)ag.994 

S(S)*1.IM 

S(«)«f.9A2 

t(7)*f.m 

S(t)«f.7IJ 

3(9)Pf «543 

0(t)«117. 

0(2)«3Sf. 

0<3>M87. 
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EF>'  SOURCE  STATEMENT 


IE'j(S) 


0<*>s \?2. 

0 ( 6 ) *1 , 

0(7)«i53^. 

0<fl)*t765. 

T(O)*9"0-». 

rpPs'i .  ,• 

Gpf’sS.  * 

DC  7  -si, 9 

"SpH«?.#c.1S(RAOO«0(J)*<!IBM) 

SnPh«7.*S!N(RA0D»0(J)#SIPH) 

GPP«Co°^C(J)*CSPH 

7  rPP«FoP*S t J)*SNPH 

8  R£TURm 
EK'O 


&9 


o ooooo 


SU811 


EFN  SOURCE  STATEMENT 


IFN(S) 


THIS  FUNCTION  EVALUATES  THE  WEIGHTED  SUM  OF  A  SERIES  OF 
BESSEL  FUNCTIONS.  H  IS  USED  TO  CALCULATE  THE  SCATTERING 
FROM  A  CYLINDER. 

COMPLEX  '’UNCTION  3ESF ( AKA » XCB •  VSR ) 

COMPLEX  SjM 

OATA  Pl.EE/3.i41592S5.2.7lS281P3/ 

CB*XC9 

SUH«(-2.2,F.) 

IFICB  .LT.  -.99996)  GO  tq  6 
SB»XSR 

I*<ABSCC9*SB)#8.*18.*AKA»1.3>/2. 

FM«I«2  . 

CB2«S0RT((l.«CB}/2.) 

V«2,*AKA*CB2 

VI«2./V 

XI*3./V 

IF ( V  .LT.  3.)  GO  TO  1 

PMlaV-.765398l6-XI»(.04166397*X!*(.00m«54-XI«{ 
..f|i!262573-Xl*(.BP054125*XI»(.|BB29333-XI«. 08813558))))) 
FO*.79708456-Xl*(.0f0f8077«XI«(.00552740*XI*( 

. .00009512-XI*( .00137237-XI«( .88072P09-XI*. 00214476 ) ) ) ) ) 
BJ«FO«COS(PHI)/SQRT(V) 

GO  TO  2 

1  XI«V*V/9. 

BJ»l.-XI*(2.2499997-Xl*(1.2«54204-XI»{.3i03844-Xl#( 

. ,0444479-XI»{ ,0039444-XI*{ .88821) ) ) ) ) ) 

2  R T ■-S0«S I N ( AKA*SB ) / C AKA* ( 1 . *CB ) ) *2 . *C82«CB2«8 J 
SER»f . 

FN«FM*10. 

E J«  < ( 1 . -1 . /FN) •• { FN- . J  > ) •EC*V/2 , /FN 
Oj»l. 

FH*FN-1l 

3  fj»-Cj*fn»vi«oj 

.  AM* ABS <  F J ) ♦ A8S ( 0 J ) 
ej*oj/am 
OjaFJ/AH 
FN«FN-1. 

IF(FN  .CT.  FM-.5)  GO  TO  3 
B«ATAn2(S0,C8) 

Sl»f In( (Fn*2. )«8/2. > 

CiaC0S((FN*2.)«B/2.) 

S2aflaCB-CiaS8 

C2«CiaC8*Sl*SB 

4  YlaFN 
2IaFN*2. 

SERaSER*EJ*(C2/YI-Cl/II ) 

IF(FH  .LT.  2.)  GO  TO  5 

ClaCS 

51*82 

TEMpac2aCB«S2«SB 

Sj«82*CB-r2»9B 

C2«TEMP 


8U  -  *FN  SOURCE  «TiTEMFl»T  -  ^(S, 

FJ*-E.j*F\'*vI»OJ 
EJ*0J 
Oj«FJ 
FMF'J-l. 

Fj*-Ej*F^*v  I  «CJ 

AW«aPs(0J)*A9S(FJ) 

FN*FN-i. 

SFR*Srff/4M 

Ej*0J/A* 

0J*FJ/4H 
CO  TO  4 

5  Aj*-Ej*FM#v|*OJ 
0J*0J*6j/4j 
SE«*Sfc#9j/aj 

R|«PI-2.*CB2*SE» 

C|»-Pl#CP2*0J 

SUH«CHPl*(RI,cI> 

6  SESF*SUH 
*£TU*W 
END 


4*4/ 
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«  r  f3W 


S'UBl? 


EFN  source  STATEMENT 


IFn<S> 


c 

c 

c 

c 

c 

c 


this  IS  The  ^ihc  function,  IT  IS  DEFINED 
X  DIVIDES  RY  X.  SINC  OF  ZERO  IS  TAKEN  TO 

FUNCTION  SINC(X) 

XX«A8S<X> 

IFCXX  .Lt.  ,fgRl22R703)  *X«,0§01 
SINC*?IN(XX)/XX 

return 

end 


AS  TWE  si Nr  or 
RE  ONE, 


SUB13 


EFi«  SOURCE  STATEMENT 


IF'KS) 


block  data 

COMPLEX  2.<H,2jP,2jPC<2).2jMC<2> 

COMMON  /AB/  2JM, 2jP,2jPC»2JMC 

COMMON  /VAR/  SM,SNCUT,SNCUD.SNCUC(?)»VPC(2)»VMC{?) 
COMMON  /$»JR/3UHM*(4>,P1,RA00 
COMMON  /ANT/0UM(43),AS(2> 
oata  A?/l. 0,1.0/ 

OaTa  SM/3S7,/ 

DATA  2JM*2JP»2JPC,2JMC  /*•<».. ?.>/ 

OATA  PI. RA00/3. 14159263.. 017413292/ 

END 


APPENDIX  B 


DYNAMIC  EMULATION 
PROGRAM  DYNM  LISTING 


The  ILSLOC  program  calculates  the  CDI  at  each  point  in 
space;  this  CDI  includes  the  Doppler  effects  from  the  velocity 
of  the  aircraft.  In  the  simulation,  the  receiver  system  is 
assumed  to  generate  the  CnI  value  instantaneously.  In  the  real 
case,  the  inertia  of  the  electrical  and  mechanical  ortions  of 
the  system  limit  the  rate  of  change  of  the  CDI.  Thus  the  real 
observed  CDI  appears  to  have  been  low-pass  filtered  from  the 
instantaneous  CDI. 

The  program  DYNM  takes  the  output  tape  generated  by  program 
ILSLOC  and  converts  it  to  observed  CDI  by  simulating  the  effect 
of  a  low-pass  filter.  The  variable  TAU  is  the  time  constant  of 
the  effective  filter.* 


Note:  When  a  flight  path  has  been  segmented,  the  low-pass 
filter  will  operate  continuously  over  the  entire 
flight  path. 
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* IpFTC  MAIN 

C  THIS  PROGRAM  SIMULATES  THE  EFFECT  OF  THE  MECHANICAL  AND  ELECTRICAL 
.C  INERTIA  OF  THE  ILS  RECEIVER  ON  THE  CDI.  THIS  EFFECT  IS  EQUIVALENT 
C  TO  A  SIMPLE  R-C  LOW  PASS  FILTER.  THE  VARIABLE  TAU  IS  THE  TIME 
C  CONSTANT  OF  THE  EFFECTIVE  FILTER.  A  TYPICAL  VALUE  IS  .4  SECONDS. 

C  THE  INPUT  TAPE  IS  ON  UNIT  11*  THE  OUTPUT  ON  UNIT  12. 

C 

C 

DIMENSION  XY ( 10 ) ,DEF { 501 ) .MEMO (14) 

LOGICAL  FOF 
DATA  ILBL/4HDYNM/ 

DATA  TAU/0.4/ 

IF(EOF(im  GO  TO  4 

1  it*o 

DELC=C . 

2  READ! 11* 1000 )  MEMO.XY, ID.NC. ICF 
WRITE( 6*1003)  MEMO.XY* ID.NC.ICF 
DEFK=ABS(XY(9)/XY(5)/TAU) 

IR=IFIX(XY(10)+.l) 

READ( 11*1001)  ( DEF ( I) » 1*1 » IR ) 

I F ( I T  .EO.  0)  CEF2*DEF ( 1 ) 

I  T*1 

DO  3  1*1, IR 

CEF2=CEF2+DELC 

DELC* ( DEF( II -CEF2 I *DEFK 

3  DEF ( I ) =CEF2 
MEMO( 13 ) *ILBL 

WRITE! 12,1000)  MEMO.XY. ID.NC. ICF 
WRITE! 12,1001)  (DEF! I) ,I-1,IR) 

IF( I D  .GT.  13)  GO  TO  1 
IF(  ID  .EQ.  0)  GO  TO  1 
GO  TO  2 

4  REWIND  11 
END  FILE  12 
REWIND  12 
CALL  Ex*.T 

1000  FORMAT (13 A6,A2,/,1X,7F  18.9, /.3F18.9, 110, 10X, 2110) 

1001  FORMAT (7E15,8) 

1003  FORMAT ( IX. 13A6  »A2» /♦ IX  »7F18.9,/»3F18.9» 1 10, 1 OX .2 1 10 ) 

STOP _ _ _ 

END 


. . . .  J.  Ji  wa»  -,.  riffcjitrtttliliMfrjMgili  “f- tv^WMlTff 


This  program  has  been  written  to  generate  graphs  of  the 
static  and  dynamic  CDI's.  It  was  written  on  the  IBM  7094  using 
the  CALCOMP  plotting  subroutines. 

The  first  input  card  has  the  following  format: 


Col. 

Symbol 

Usage 

1-2 

NL 

Number  of  lines  per  graph 

3-4 

NGRFS 

Number  of  graphs 

5-7 

NTAPE  (1) 

Input  logical  unit  no.  for 
line 

first 

8-10 

NTAPE  (2) 

Input  logical  unit  no.  for 
line 

second 

11-13 

NTAPE  (3) 

Input  logical  unit  no.  for 
line 

third 

NL  pe-  mits  the  overlaying  of  two  or  more  CDI  or  signal 
strength  graphs  for  comparison  purposes.  The  scaling  will  be 
set  by  the  first  graph,  and  the  successive  overlays  will  be 
plotted  to  the  same  scale.  A  maximum  of  three  lines  per  graph 
will  be  allowed. 

NGRFS  set. 5  the  maximum  number  of  graphs  to  be  drawn.  Each 
graph  will  have  the  same  number  of  overlays. 

NTAPE  (i)  gives  the  logical  unit  number  used  for  the  input 
of  the  ith  line  on  each  graph.  If  the  value  of  NTAPE  is  negative 
then  its  absolute  value  will  be  used  as  its  logical  unit  number 
and  the  tape  will  be  rewound  before  input. 

The  second  input  card  defines  the  scaling  used  for  the  graph 
(or  graphs)  described  above.  It  has  the  following  format: 


m 


t'±f_n  ■wi^f«Thc 


i^"j=U5S3^ 


ft 


Col .  Symbol 

I- 10  XSC 

II- 20  DELX 

21-30  YMAX 

;  ■  .)  YMIN 

si -SO  DELY 


Usage 

Horizontal  scale  in  ft. /in.  or  deg. /in 

Tick  mark  spacing  in  ft.  or  deg. 

Maximum  y-value  on  vertical  scale 

Minimum  y-value  on  vertical  scale 

Tick  mark  spacing  on  vertical  spacing 
in  microamps  for  CDI  or  relative  units 


The  horizontal  axis  is  drawn  in  either  feet  or  degrees  per 
inch  as  specified  by  XSC.  The  tick  mark  spacing  along  the  axis 
is  determined  by  DELX.  The  length  of  the  axis  will  be  adjusted 
to  the  shortest  length  with  an  integral  number  of  tick  marks 
that  will  cover  the  domain  required  by  the  input  data.  When  a 
flight  path  has  been  segmented  it  is  treated  as  a  single  line 
on  the  graph. 


YMAX,  YMIN  define  the  range  of  the  plotted  variable:  CDI 
or  relative  signal  strength.  The  Y-axis  has  a  fixed  length  of 
seven  inches.  If  DELY  does  not  integrally  divide  the  range,  DELY 
will  be  adjusted  to  yield  an  integer.  When  the  range  (YMAX-YMIN) 
is  zero,  the  program  will  automatically  scale  the  range  to  the 
largest  scale  that  will  include  the  data  in  the  length  of  the 
axis . 


When  multiple  graphs  are  plotted,  each  graph  is  scaled  in¬ 
dependently. 

After  all  NGRFS  graphs  have  been  drawn,  the  program  will  loop 
back  to  the  beginning  and  attempt  to  read  in  a  new  NL  card.  This 
allows  many  graphs  to  be  drawn.  If  the  user  wishes  to  replot 
data  using  different  scales  or  overlaid  with  different  sets  of 
data,  he  may  use  the  negative  NTAPE  to  rewind  the  input  tape. 

The  program  will  terminate  after  reaching  an  end-of-file 
on  the  card  input  unit. 

The  vertical  scale  on  the  graph  is  always  labeled  "micro¬ 
amperes".  This  is  valid  only  for  CDI  graphs.  All  others  are 
in  relative  units  and  this  labeling  should  be  ignored. 
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04/l< 

MAIN  -  EFN  SOURCE  STATEMENT  -  IFNI  Si  - 

COMMCN/T6ST/XMIN,OXR,NTOTtNP 
LOGICAL  EOF 
DIMENSION  I8UFI1000) 

01 MENSICN  NTAPEI3) .MEMO! 14) , Ml  14) 

EQUIVALENCE  (M( 1) .MEMO! 1) ) 

CCMMCN  /PDF/  CFI2000 )»XLFN, NSTEPS, IDEF, I DENT *0X1  IQ' »NPTS( 10) 

COMMON  /PRINT/  NL ,XSC ,DELX ,YMAX, VMIN, OELY, ICF 
CALL  PLOTS! IBUF, 1000) 

CALL  FL0T(P.0,-12.,-3) 

CALL  FACTOR  10.4) 

IL8L*l 

60  CONTINUE 

I  FI  EOF! 5))  GO  TO  55 
RFADI5, 100)  NL.NGRFS.NTAPE 
WRITEI6* 100)  NL,NGRFS#NTAPE 
I  FI NGRFS.LE. C)  NGRF  S*3 
100  FORMAT (212*313) 

DO  401  1*1. NL 

IFINTAPEID.GE.O)  GO  TO  401 
NTAP6II)*-NTAPEII) 

NU=NTAPE( I ) 

REWIND  NU 
401  CONTINUE 

READ! 5 . 1G1 )  XSC .DFLX.YMAX.YMIN.OELY 
WRITEI6. 101)  XSC .DELXf YMAX. YMIN.OELV 
101  FORMATC8F10.0) 

TEMP*AMIN1 I Y MIN. YMAX) 

YMAX*AMAX1I YMIN, YMAX) 

YMIN-TEMP 

TEMP-YMAX-YMIN 

IFITEMP  .NE.  0.)  DEL Y* TEMP/!  FLOAT!  IFI XI TEMP/QELY*.  5)  ) ) 

NPLT  *  1 
NP  *  l 
I  =  l 
Nl  *  1 
NTOT  *  0 

10  NU  =  NTAPEINP) 

IFIEOF(NU) )  GO  TO  50 

READINU. 500)  M. XO.OXR.XY.ID.IOEF.IDENT.ICF 
IFIICF  .NE.  0)  ICF* 1 

WRITEI6.600)  MEMO.XO vDXR.XY. 10. IDEF, IOENT, ICF 

IFIILBL  .NE.  ii  GO  TO  70 

IL8L*0 

CALL  SYMBOL 10..0. ..14.MEM0.90..80) 

CALL  PLOTI3..0..-3) 

70  CONTINUE 

IR  *IFIXI  XY+.l) 

NTOT  *  NTOT  ♦  IR 
IFII.E0.1)  XMIN  *  XO 
500  FORMAT I 13A6.A2././.3F 18. 9.4110) 

600  FORMAT I2X. 13A6.A2./. 3F 18.9, AI 10) 

501  FORMAT  I 7E15.8) 

502  FORMAT I 1X.7E15. 8) 

READ (NU, 501 ) I DF! J ) » J*Nl,NTb • ) 

WRI TE (6,502)  I DFI J) , J*N1 ,NTCT ) 


70- 


4 

1 


*AIN 


-  EFN  SOURCE  STATEMENT  - 


I FN  ( S  )  - 


•4/1 


10001  XMIM.IR.Nl.NTOT.NP,! 
1000  FORMAT (F 10. 0,5110) 

N°TS( I I  =  IR 
Dxm  =  dxr 

IF(  IC  .GT.  13  )  GO  TO  40 
IF(I0  .SQ.  C)  GO  TO  40 
M  *  N1  ♦  JR 
1  =  1  +  1 
CO  TO  10 
11  NL  =  NP 
40  CONTINUE 
NSTcps  =  I 

IF(NP.GT.l)  GO  TO  41 
CALL  GRAPH2I0I 
GO  TO  42 

*1  CALL  GRAPF2C 1  ) 

42  continue 

N 1  =  1 
I  =  1 
NT  IT  *  0 

I F( NP.EC.NL )  GO  TO  45 
NP  =  NP  ♦  1 
GO  TO  10 
45  NP  =  1 

CALL  PLCT(XLEN  +  7.,-)2.,-3) 

NPLT  =  NPLT  ♦  l 
IL3L*  1 

IF  I NPLT .GT.NGRFS  )  GO  TO  60 
GO  TO  10 
50  CONTINUF 

IFINTOT cGT.O )  GO  TO  11 
CALL  PLOT  ( XLEN+7. ,-12. »— 3 ) 

GO  TO  60 
55  CONTINUE 

CALL  PLCTIO. ,0. ,999) 

DO  400  1=1, NL 
NU=NTAPE II) 

4C0  REWIND  NU 
STOP 
PNO 
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SUB1  -  EFN  SOURCE  STATEMENT  -  IFNIS)  - 

SUBROUTINE  GRAPH21 ITL ) 

DIMENSION  XLA8C4) 

C0MM0N/T6ST/X0»DFLTAX»MDELTA ,NP 

OAT  A  XLA0/24HOISTANCE,FT.  DEGREES  / 

DIMENSION  TYPEI3I 
OIMENSION  X( 3) ,NC! 3) 

COMMON  /PDF/  CF! 2000 ),  XLEN, NSTEPS* IDEF , IDENT»DX( 10) »NPTS!  10) 
COMMON  /PRINT/  NL ,XSC ,DELX ,YMAX, YMIN, OELY, ICF 
CATA  X  /-5.» 5.,  5./ 

OATA  NC  / 1 » 5»4/ 

IFC ITL  .NE.  C)  GO  TO  I 
ELX*D£LX 

Ic( DELTAX.LT. 0. )  ELX  *  -ABS(OELX) 

RANGE*0. 

DO  11  1*1, NSTEPS 

RANGE*RANGE*FLOAT INPTS II)  )*CXl I ) 

TI X* I FIXC RANGE /EL X*. 9) 

XLFN  *  A8S(ELX/XSC*TIX) 

IFIXLEN  .GT.  40.)  GO  TO  9 
1FIXLEN  .GT.  5.1  GO  TO  6 
9  XSC*AeS IRANGE/20. ) 

XLEM*AbS(ELX/XSC*TIX) 

WRITE! 6, 8)  XSC 

.o?^.4,Ti25H  Axls  0UT  0F  RANGE  SCALE* ,E  12.  5,  8H  FT./IN.  /) 
continue 

XMAX*TIX»ELX*X0 
XMIN  =  APIN1IX0,XMAX) 

XMAX  *  AMAX1 t XO ,XMAX) 

ND  *  2 
PWR  *  0. 

CALL  PLOT (0. ,1. 5,-31 

AMIN*YMIN 

AMAX*YMAX 

rAiVivvcstA  M*ALL^SCLAX(  7*  »DF*NOELTA,AHAX, AMIN, OELY, NO, PWR) 

YSCL*4OELN<0*,0*,  AX,AMIN,0ELV,7*,l2HMICPOAHPERES,l^’N0,P-R,0ELN, 

IXLAB*2*ICF+1 

IXSC*-1 

IF! ABS1CLX1  .LT.  10.)  IXSC*1 

Cn=!\11X*S3<0*,0**XHA*,XMIN,E*’X,”XLEN,XLAB*IXLAB)  *12  fIXSCfO. 

•  t  UcLN ) 

XSC  *  OPLN 

XT  =  XLEN/2.  -  2. 

I  FI  AM IN*AMAX .GT .0.1  GO  TO  2 
Ip!  AMIN  .EQ.  0.)  GO  TO  2 

ZERO* !0. -AMIN/ 10.**PWR)/YSC 
CALL  PLOT  !0«, ZERO, 3) 

CALL  PLOT ( XLEN, ZERO ,2) 

CONTINUE 

CONTINUE 

XI*0. 

IF!OELTAX  .LT.  0.)  XI»XMAX-XMIN 
J*l 

DO  5  1*1, NSTEPS 
OELTAX  *  OX! I) 
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SUPl  -  EFN  SOURCE  STATc»,ENT  »  IFN(S) 

NX=NPTSCI) 

IF(I  .LT.  NSTf PS )  NX=NX* 1 

ym=apin/io.-*apwr 

CALL  XCLINF(XI»0£LTAX*DF(J| ,NX,0. ,XSC,YM,YSC,NC (NP) ) 
•l*J  *-NPTS  ( I  ) 

XI=XUCX(I)*FLCAT(NPTS(I>) 

5  CONTINUE 
RETURN 
END 


^ 
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SUB3  -  EFN  SOURCE  STATEMENT  -  IFN(S)  - 

SUBROUTINE  XCLINECXI ,DX,Y,N,XF,DELX,YM, DELY,NC ) 

DIMENSION  Y<1)*IPEN<4) 

PEAL  L (4*4 )*LL(4i 
DATA  I  PEN/2 # 3*2 *3/ 

DATA  L/*3» • 1*  *3»* 1 *• 5*3*# 05* • 3 » 3** lt*l»*05*«lt* 05/ 

X  =  XI 

2  IC  *  NC  -  l 

XPl  *  (X-XMJ/DEIX 
YP1*(Y(1)-YN)/DELY 
CALL  PLOT ( XPl* YP1 *3) 

IF< IC.LE.OJ  GO  TO  10CO 
IFIIC.GT.4)  IC  *  4 
K=1 
1*2 

X  *  X  ♦  DX 
XP2  *  ( X— XM) /DELX 
YP2*  (Y(2)«YM)/DELY 
1  LL( K) »L(K *  I C ) 

10  DIFFX*XP2-XP1 
0JFFY*YP2-YP1 

DIS*S0RT|0IFFX*0IFFX40IFFY*DIFFYI 

IF(DIS.GT.LL(K) I GO  TO  100 

CALL  PLGTCXP2*YP2»IP£N(K)I 

XPl*XP2 

YP  1*YP2 

1*141 

I  FU.GT.N3  RETURN 

X  *  X  ♦  DX 

XP2  *  (X-XMI/OELX 

YP2*IY! I l—YMI/OELY 

LLIK)*LLIK)-DIS 

GO  TO  10 

IOC  RATIO*DI S/LL(K) 

XP1-XPUDIFFX/RATI0 
YP1*YP1*DIFFY/RATI0 
CALL  PLOTIXPlfYPlflPEN(KI) 

K»K>1 

I F  € K.EC.5) K*1 
GO  TO  1 

1000  DO  50  I*2*N 
X  *  X  ♦  OX 
XPl  *  ( X— XM) /DELX 
YPl* (Y( I J-YM J/DELY 
5C  CALL  PLCTCXPl,YPi,2l 
RETURN 
END 
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SCLX  -  FF*  SOURCE  STATEMENT  -  IFN(S)  - 

S'liJc OUT! NP  SCLAX<AINCH,VAa,N,VMAX,VMlN,DEL7A,ND, EXP) 

OIMSNSIGN  VAR(l) 

AXLEN  =  AIKCH 
VM  AX  =  VAM1J 

VII  N  *  V  A  ?» ( 1 )  *, 

00  40  1*2, N 

VMAX  *  AMAX1(VMAX, VAR{ 1 ) ) 

40  VMIN  =  AWINllVMINtVAAd)  I 
MO  *  C 

ME  =  0 
M  =  2 

TOTAL  =  VMAX  -  VMIN 

DETERMINE  EXPONENT  ANC  INCREMENT/ INCH 
VM  =  fiNAXl (ABS<VMAX),ABSIVMIN) ) 

IF(VMAX*VM!M  6,5,7 
7  VAV  =  ARSC VMAX+VMINJ/2. 

DELTA  *  T'JT AL/AXLFM 

IPC  TOTAL. GT. 0. .AMD. TCTAL/VM.LT . . 751  GC  TO  4 
IFCVMAX.E J.VM)  V«IN=0. 

IF(VMIN.PJ.-VM)  vmax*o. 

GO  TO  5 

6  AXLEN  *  AXL6N*VM/T0TAL 
5  DELTA  *  V*/ AXLEN 
WAV  *  VM/7. 

TEST  FOR  VAV  BETWEEN  .01  AND  1000. 

4  IF(VAV.Lt. 1. £— 1 1 )  GC  TO  21 
I F C VAV  -  .01)  3,10,1 

41  IFCVAV  -  1.)  3,10,10 

1  IF ( VAV  -  10DC.)  10,2,2 

VAV  GE  1CCC. 

2  IFtNE.EC.D)  VAV  *  VM 
VAV  =  VAV/100U. 

N£  *  NE  -  3 
GO  TO  l 

VAV  LT  1. 

3  VAV  =  VAV*10G0. 

NE  *  NE  ♦  3 

GO  TO  41 

DETERMINE  DECIMAL  PLACES  IN  DELTA 
1C  IFIDELTA.LT. V-/1.E4)  GO  TO  21 
DELTA  *  0ELTA*10.**NF 

11  IFIOELTA  -  1.)  12,19,13 

12  OELTA  *  CtLTA*10. 

NO  *  NO  «■  l 

GO  TO  11 

13  IFtOELTA  -  10.)  15,8,14 

14  DELTA  =  DELTA/iC. 

NO  *  NO  -  1 

GO  TC  13 

DELTA  NOW  BETWEEN  1  AND  10 

15  IF(OELTA  -  5.)  16,17,17 

16  IFIOPLTA  -  2.)  19,16,18 

17  DELTA  a  5./lC.**(ND*NE) 

GO  TC  20 
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SCLX  -  EFN  SOURCE  STATEMENT  -  IFNI  SI  - 

Id  DELTA  *  2./10.**CND*NE) 

M  *  5 
GO  TO  20 
8  NO  «  NO  -  l 

19  OELTA  »  l./10.**IND*NE) 

C  RESET  VN IN  t  FIRSTVi  FOR  AXIS 

20  AK  *  VMIN/OELTA  ♦  .01 
K  *  (  I F IX (  AK  I/M  t  *M 
IFCVMIN.LT.O.)  K*K-H 
VMIN  *  0ELTA*FL3ATIK1 

NOIV  *  (VNAX  -  VH IN l/OELTA  ♦  .9 

IF  I  FLOAT! NOI VI .GT .AI NCH*2. 1  OELT A*0ELTA*AMAX1 (2. »  FLOAT! N)/2. 1 
IF( NO.LEoO I  NO  *  -1 

21  EXP  =  NF 

WRITEC6.1002)  VMAX,VPI N.OELTA.ND.NE 
RETURN 

1002  FD»MATC1H0.3E13.3»3I7//1 
END 
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SJdAOUTlNF  AXIS3CX0,Y0,VMAX.V«IN,nELV,AINCH,8CD,NCR,NDt:C,PW&,VSC) 
FACTOR  =  m.**PWP. 

AM! N  =  VMIN*FACTUR 
AM4X  =  vmax*factcr 
OsLX  =  ABS  C  CELV  MFACTOP. 

0IMFNS1CN  BCOC1) 

H7  =  .15 

*1-0. 

12=0. 

Ml  *  0. 

NEXP  *  0 
NCH=IABSCNCR) 

IF(PWR.N£.3. )  NCXP  =  6 
CI\CH=i»eStAlNCH) 

I F  (  f  VMAX— V**I  K }  /  A MAXI  I VMAX, ~VWI N1  .  LT  .  1 . £-61  GO  TC  50 
IF({AMAX~AMIN)/(DFLX«-1.F-8).GT.>3.*CINCH)  OFLX  =  •  t  AMAX-AMIN J/CI  NCH 
IF ( DELX.GT  .A  “AX- AM  IN  )  DELX  *  AMAX  -  AMIN 
IF ( NCR.LT. 01  W3  -  1. 

VIM* { AMAX-AMIN ) /0ELX*1. 9 
4NC*C!  NCM/FLGAT  1 ) 

IFIAINCH.LT. C.JGO  TO  5 
W2*l. 

GO  TO  10 
5  Wl*  1. 

13  CALL  FLCT CX3 » Y0»  3 1 
VSC  *  OELX/FACTGR/ANC 
AN'JM*A«IN-OELX 
X*3. 
r=o. 

XM=0. 

OFF  *  .05 
on  AO  1*1, Mjp 
4NJM*ANUH+DeLX 
11*3 

25  IF C  ABSC  AKUM) /10.**IJ.LT.1.1G0  TO  20 
11*11+1 
GO  TC  25 

20  I  F( ANUP.LT. C. 111*11+1 

ifcabscanuhi.lt. 1.)  11*11+1 
iMORp*Nrec«-i 
1 1  =  1 I+IMORE 

IF1  IFIXCWlMI.CQ.il  HT  *  AMIN1CHT  ,ANC  /FLGATC I  M2 1 1 

ML  *  AMAXl C.12,1.2*HT1 

CCmTEP  =  FLCATCIIMHT/Cl.  +  WiJ 

XC  *  X  -  CFNTCR  -  W2*.15 

IFCXC.LT.XP)  XH  *  XC 

I F C  W2*W3.GT.C.l  XC  *  .15 

IFCABSCXCl.GT.ABSCXHli  XH  =  XC 

YC  *  Y  -  Ml*  CHT  +  .15  -  M3+CHT+.31)  -  W2+0FF 

CALL  PLOT! >C+X« Y0+Y,2' 

CALL  PLOTtX3+X+.l*W2»YO+Y+.l*Wl»21 

CALL  PLOT! XO+X-. 1**2 ,YO+Y-. 1+W1 ,2 1 

CAL L  NUHBFq (  xO+XC ,  YO+YC , HT, ANUN, 0.  ,NO£C  1 

CALL  FLCTCXO+X,VO+Y, 31 

X=X+4NC*hl 
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AST  =  (CINCH  -  FLOAT(NCH*NEXP)*HLJ /2. 

IF*'  ^"'•1#)  XM  s  -XH 

*  --  ;  (XO  *  BST )  ♦  W2*(XO  ♦  XF  -  OFF  +  W3*(2.*OFF+HL ) I 

'  •  »,•>  (YO  *■  VC  -  1.5*HL  +  W3*(HT  *  2.*HLM  ♦  W2*( YO+BST > 

.  v.  J  N  ...  jPL(  XXC#YYC»HLt  BCD»9C.*W2,NCH) 

\,..EQ.  O.  )  RETURN 

CALL  SYM8CL(999.,999.,HL,5H  *  1C,9C.*W2,5) 

X  =  999.  ♦  (XXC-.6A*HL-999.)*W2 

Y  =  999.  *  <YYC*.66*HL-999.)*W1 

CALL  NUMBER ( X, Yt . 75*HL» PWR » 90.*W2»- II 
RETURN 

5C  VSC  =  <  VMAX-VMIN«-l.E-6/FACYCR)/CINCH 
WRITE (6, 1000) 

10UC  FORMAT ( IH9 »27HINSUFFICI ENT  RANGE  FOR  AXIS  ‘ 

RETURN 

FND 
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